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TECHNICAL MEMORANDUM X- 64727 


A COMPARISON OF CMG STEERING LAWS FOR HIGH ENERGY 
ASTRONOMY OBSERVATORIES (HEAOs) 

SECTION I. INTRODUCTION 


During the past year and one-half, the Preliminary Design Office, 
Program Development at Marshall Space Flight Center (MSFC) has conducted 
rather extensive Phase A studies of spacecraft designed to detect and observe 
high energy radiation sources. These spacecraft have been designated as high 
energy astronomy observatories (HEAOs) . Missions A and B are supposed 
to scan the entire celestial sphere over an extended time period and, then, 
point to selected radiation sources [ 1] . Mission C has been designated only to 
point to selected targets. All HEAO configurations utilize solar panels to 
receive power and, hence, must be solar oriented within certain power and 
thermal constraints. HEAO-C, however, has more demanding pointing specifi- 
cations than HEAO-A or -B. Due to the limitations of fuel weight and inherent 
limitations in pointing performance of an all thruster reaction control system 
(RCS), control moment gyros (CMGs) have been baselined for HEAO-C. 

An RCS will be used for momentum management of the CMGs, but electro- 
magnets offer more growth potential and allow continuous CMG momentum 
dump without interrupting vehicle pointing (Appendix A). This report contains 
a description of the CMG system which has been selected during preliminary 
studies for HEAO. More specifically, the orientation of the CMGs relative 
to the HEAO reference axes have been selected to provide a near spherical 
momentum envelope with all CMGs operational and also provide complete 
vehicle control, even with one CMG failed. 

Four single gimbal (SG) CMGs are arranged in a skewed configuration 
about the sun pointing vehicle axis. To provide effective vehicle control 
torques, the CMGs must be gimbaled in response to an attitude error signal. 

As will be shown, the gimbal commands are not necessarily unique but depend 
upon the assumptions made to obtain a solution to the CMG torque equations. 
Whatever the solution, it is referred to as the CMG "steering law.” Several 
condidate steering laws are derived and evaluated according to their effective- 
ness in producing the control torque required by the attitude error signal. 



The candidate steering laws are the constant gain, MSFC maximum 
contribution, pseudo inverse, Bendix three gimbal inverse, General Electric 
(G. E.) transpose with torque feedback, and the Teledyne Brown Engineering 
Company (BECO) H-distribution, Each steering law was utilized in conjunc- 
tion with the Euler equations for HEAO (Appendix B) with the four skewed 
SG CMGs in a digital simulation (Appendix C) to obtain the corresponding 
vehicle pointing performance. Each steering law was then evaluated on the 
basis of complexity in implementation, accuracy ol pointing performance, 
avoidance of mathematical singularities, possible CJVLG gimbal angle positions 
which prevent the desired torque from being produced (gyro hang-up) , adap- 
tion to a CMG failure, and performance after a failure. Based on accumulated 
study results, the pseudo inverse CMG steering law is recommended for 
HEAO. 


SECTION II. CMG CONFIGURATION SELECTION 
A. Introduction to CMG Control 

The path of an orbiting vehicle is predetermined by its orbital param- 
eters, essentially that of a free-falling ballistic trajectory in its orbit. The 
basic objective of the spacecraft attitude control system is not to change its 
orbital path but to maintain a prescribed attitude (orientation) as a function 
of time relative to inertial space irrespective of the flight path. The principle 
of conservation of angular momentum led to the "momentum exchange" idea, 
whereby in the absence of externally applied torques, if one part of a closed 
system increased its momentum by a specified amount, the remainder of the 
system lost an equal amount of momentum. An example is a flywheel sup- 
ported by a frame hard-mounted to the spacecraft with a momentum given by 


where 1^. is the ilywheel inertia, is the wheel angular velocity, and 

is the Ilywheel momentum. By decreasing the flywheel speed, a torque is 

generated about the flywheel spin axis which counter rotates the spacecraft. 

The angular momentum stored in the Ilywheel decreases in proportion to the 

change in flywheel velocity while the spacecraft momentum H increases 

v ’ 
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Since the flywheel support frame is hard-mounted to the spacecraft, the torque 
applied to the spacecraft is given by 



After the flywheel momentum has been decreased by a prescribed amount, 

AWf , the total momentum of the flywheel and spacecraft must remain 

constant assuming that no external torques act on the spacecraft. The change 
in spacecraft velocity, AW^ , due to an arbitrary change in flywheel momen- 
tum is 


AW 

v 


AW 


I 

v 


(3) 


where I is the vehicle moment of inertia about the axis aligned with the fly- 
wheel spin vector. 

Consider a constant-speed flywheel that is mounted on a gimbal rela- 
tive to the spacecraft. By rotating about the gimbal axis, the spin axis of the 
flywheel and the direction of its momentum are changed relative to the vehicle. 
Although the flywheel momentum remains constant, a gimbal rate a produces 
the torque 


T = -(a x H ) ' (4 ) 

which is perpendicular to both the gimbal axis and flywheel momentum vector. 

After the flywheel momentum has been rotated by a prescribed small 
amount, Aq , the corresponding change in the spacecraft velocity required 
to conserve angular momentum can be approximated by 

Aq I W 

*W*-— I • (5) 

v 
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In the first example, the flywheel momentum value was varied, 
whereas, in the second example the orientation of the flywheel momentum 
vector relative to the vehicle was varied to obtain a toi'que. For either case, 
the vehicle is reoriented relative to inertial space so that momentum is con- 
served assuming no external disturbances. The first method is called reaction 
wheel control, while the second is called gyro control. Generally, the gyro 
rotor has constant speed of rotation and is referred to as a control moment 
gyro. If only one gimbal is used to reorient the momentum, the CMG is 
referred to as a SG CMG, whereas, if two gimbals are used, the CMG is 
referred to as a double gimbal (DG) CMG. Since the vehicle requires three 
degrees of fieedom to maintain a prescribed orientation relative to inertial 
space, the attitude control system must have actuators that provide effective 
control torque about three independent axes. Hence, a momentum exchange 
system must provide at least three degrees of freedom for control purposes. 
Conceivably, one variable speed reaction wheel mounted on two gimbals could 
provide three-axis control for the spacecraft. If SG CMGs are utilized for 
momentum exchange, then at least three units are required for three-axis 
spacecraft control. In addition, the three units must be mounted relative to 
each other such that three independent degrees of freedom are obtained. 


B. Single Gimbal CMG 

The characteristics oi a SG CMG are illustrated in Figure 1. The 
flywheel turns at a constant speed producing momentum directed along the 
V axis - The flywheel momentum is rotated in the Y^-Z plane by gimbaling 

about the X-axis, thereby producing a torque along the Z^-axis. The torque 

produced obeys the vector cross product law, equation (4). Therefore, at any 
instant of time, the torque produced must be perpendicular to both the gimbal 
axis and the momentum axis. For example, assume that it is desired to 

generate a vehicle torque which lies in the X -Y plane. Then the desired 

0 0 

torque cannot be generated by the SG CMG shown. For any desired torque, 
only the piojection of that torque onto the Z —axis can be generated. From 

this simplified example, it is apparent that at least three SG CMGs must be 
utilized to obtain effective three-axis vehicle control. Moreover, the three 
CMG torque axes, Z^(c - 1, 2, :i), must span a three-dimensional vector 

space. When the torque axes oi a SG CMG system are coplanar, a control 
torque perpendicular to that plane cannot be produced. Such a condition is 
referred to as gyro hang-up (ghu). 
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Figure 1. Single gimbal CMG. 

The first general problem area is the selection of a momentum exchange 
system that is appropriately sized to counteract the environmental disturbance 
torques and, in addition, satisfy all specified vehicle maneuvering require- 
ments. 


The selected system must provide the following: 

• Sufficient reliability/ redundancy over the mission duration. 

• Sufficient torque to counteract disturbances. 

• Adequate momentum storage. 

• Adequate maneuver rates. 
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• Spacecraft stability. 

• Desired response characteristics. 

• Adequate degrees of freedom. 

As general design criteria, the CMGs should provide enough momentum 
to counteract all environmental disturbances over a one-orbit period before 
desaturation is required. In addition to cyclic disturbances, gravity gradient 
torque almost always produces a secular momentum component which even- 
tually saturates the CMGs. That is, the CMG system produces all the 
momentum it can in a given direction until no more can be produced. The 
mounting of the CMGs relative to the vehicle reference axes determines the 
shape of the maximum momentum envelope within which the CMG system can 
provide momentum. In general, the momentum envelope is shaped propor- 
tional to the vehicle moment of inertia values, especially for an inertially 
oriented spacecraft that does not maneuver very often. However, when the 
spacecraft is reoriented, stored CMG momentum is transferred from one 
axis to another. Therefore, the momentum envelope should be spherical for 
spacecraft such as HEAO-C where many maneuvers are made, or for HEAO-A 
where the spacecraft spins to scan the celestial sphere. Moreover, the CMG 
mounting arrangement must permit the CMG torque vectors to span a three- 
dimensional space to obtain the degrees of freedom required to control the 
spacecraft. When four SG CMGs are used, as dictated by reliability consid- 
erations for example, the most independency between CMGs [ 2] can be 
obtained by arranging the CMGs symmetrically about a vehicle axis as shown 
in Figure 2. The four SG CMGs are shown gt a zero momentum state (null 
position) and the CMG gimbal axes subtend an angle (i (skew angle) relative 
to the body reference axis X . The skew angle can be used to shape the 

momentum envelope. The mounting arrangement shown in Figure 2 is 
referred to as four skewed CMGs and has been recommended for use on the 
HEAO spacecraft [3] . 

Based on previous study results (4, 5] , CMGs offer several advantages 
over reaction wheels especially from a power and weight viewpoint. More- 
over, based on hardware availability [6, 7] , there are several SG CMGs that 
are sized appropriately for the HEAO-C spacecraft, from both a torque and 
momentum viewpoint. For these reasons, SG CMGs have been baselined for 
HEAO-C. To provide continued operation capability when one CMG fails, at 
least four CMGs must be utilized. However, more than four may be dictated 
by reliability considerations to achieve the required tw'o-year lifetime. The 
CMG system selected by Bendix for the HEAO-A [ 6] seems to satisfy the 
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Figure 2. CMG mounting arrangement relative to HEAO 

reference axes, 

HEAO-C requirements and, for commonality between the HEAO-A and -C 
spacecrafts, it has been baselined as the HEAO-C momentum exchange system. 
Figure 3 illustrates the CMG arrangement relative to vehicle reference axes. 
Each CMG momentum vector is restricted to a plane that is skewed relative 
to the vehicle Y^-Z^ plane by the angle /3 ; the four planes form a pyramid 

whose apex is aligned with the vehicle X - axis; and each gimbal axis, X , 

r c * 
is perpendicular to its associated plane as shown. The configuration is 

symmetrically skewed about the X^ vehicle axis so that none of the gimbal 
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axes are parallel and none are parallel to a vehicle axis. As a result, each 
CMG can contribute momentum along each axis of the vehicle. If one CMG 
fails, the remaining three CMGs provide the three degrees of freedom 
required for attitude control. 

Once the CMG configuration has been selected, the second general 
problem area is closure ol the attitude control loop through the momentum 
exchange system by gimbaling the CMGs in response to the attitude error 
signals. The logic and error signals which are used to drive the CMG gimbals 
are defined as the CMG steering law. The steering law must be selected such 
that the CMG torque produced closely approximates the desired vehicle control 
torque that is needed to maintain the vehicle's specified orientation. The first 
task that must be done prior to deriving a CMG steering law is to relate the 
CMG momentum and torque to the vehicle control axis. The momentum of 
each CMG must be projected into body control axes and summed to obtain the 
total CMG system momentum. In carrying out the required operations, several 
coordinate systems must be defined. 


C. CMG Reference Systems 

For any single gimbaled CMG, a coordinate system in which the CMG 
momentum is always constant along one axis (Fig. 4) is defined as follows: 

i unit vector along the gimbal axis X 
c c 

j unit vector along the momentum axis Y 
c c 

k unit vector along the torque axis Z 
c c 

The CMG coordinate system moves as the gimbal is varied with 
respect to the spacecraft body axis. Therefore, the momentum is always 

aligned with the Y -axis and the gimbal rate vector with the X -axis The 
c c 

torque produced by the cth CMG obeys the vector cross product law and always 

is aligned with the Z -axis. In the CMG constant momentum system, the 

gimbal rate A , momentum h , and the torque h can be written in 
c C C 

vector form as follows: 


a = a l 

c c c 


( 6 ) 
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( 7 ) 


h = h j 
c c c 


and 


h = axh = a h (c x j ) = a h k 
c c ccc c ccc 


( 8 ) 


A second CMG system is defined by setting the CMG gimbal angle to 

zero or to a position which nulls out the total momentum of all CMGs. Such 

a reference, illustrated in Figure 5 by X , Y , Z , is defined as the CMG 

n n n 

null coordinate system. When the gimbal angle is zero, the CMG null system 
is identical to the CMG constant momentum system. The CMG null system is 
related to the constant momentum system by the rotation a about the gimbal 
axis which is constant in either system. The subscript n denotes the null 
coordinate system for a particular CMG. The vector matrix form the trans- 
formation between the two systems is written as: 


X - A X 
c cn n 


(9) 




Figure 4. CMG coordinate system. Figure 5. CMG null coordinate. 
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where 



The manner in which the two CMG reference systems have been defined 
permits the matrix -A ^ ^ to hold for any single-degree-of-freedom CMG. 

However, the mounting of each CMG is unique. 

Each CMG has its own null coordinate system uniquely defined relative 
to the spacecraft body axis by its mounting arrangement. For each CMG, a 
matrix transformation must be derived to relate the spacecraft reference 

axis to the CMG null coordinates. The relation may be written as 


X - A X 
n nr r 


( 10 ) 


where the subscript r denotes the body axis reference frame. The relation 
between body and CMG constant momentum systems is obtained by 


X = A X 
c cn n 


(A A ) X G X 
cn nr r cr r 


(id 


where 


c 

C 

C 

§n 

§12 

§13 

c 

C 

C 

§21 

§22 

§23 

C 

C 

C 

§31 

§32 

§33 


11 


The elements of G are obtained by matrix multiplication of A and A 

cr cn nr 

and must be derived for each CMG. The letter c would take on the number 

assigned to a specific CMG. Since the transformations in this case are 

orthogonal, the inverse is identical to the transpose, which is denoted by an 

asterisk superscript; hence, 


X = G X . (12) 

r cr c 

Use of transformation (11) yields the following equations for the cth CMG 
gimbal rate, momentum, and torque in body axes: 

“c = °c (g “ *r + g * 2 j r + gl ° 3 V 
h c = h c (g 2 i i r + g22 j r + S23 k r ) 

and 


(13) 
( 14) 


• « 

h c = « c h c (gsi i r + g32 j r + S33 k r ) • ( 15) 

The equations for total momentum and torque from m CMGs is obtained by 
summing the vector components: 

m 


H (CMG) = ^ h 

-l C 

c - 1 


h i + h i + h k 
x r y r z r 


( 16 ) 


and 


m ^ 

H (CMG) - V h 
- j . c 
c~l 


hi+hj+hk 
x r y r z r 


(17) 
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Due to environmental forces acting on an orbiting spacecraft, the CMG 
momentum vectors will deviate considerably from their null positions. For 
most orbits in which the spacecraft is inertially oriented, momentum tends to 
accumulate in some direction due to biased environmental forces. Under these 
conditions the CMG momentum becomes concentrated in this direction until no 
further momentum can be obtained from the CMG system. This condition is 
referred to as CMG saturation. To desaturate the CMGs, a torque must be 
applied to the vehicle such that the CMGs are driven back either to their null 
position or some bias level by trying to counteract the applied torque. 


D. Four-Skewed CMG Configuration 

To develop a CMG steering law, the transformations, equation ( 11) . 
must be derived lor each CMG which relates its torque and momentum to 
spacecraft reference axes. The four-skewed CMG configuration, baselined 
lor HEAO-C, is illustrated in Figure 3. Each CMG is shown at its null posi- 
tion and the geometry between the CMG null and spacecraft reference coordi- 
nates is illustrated. At the null position the momentum of CMG number 1 and 
number 3, as well as that of CMG number 2 and number 4, cancel. The 
transformations are carried out by first rotating negatively about each Y 

n 

axis by the angle /3 which aligns the transformed X ' axis with X reference 

n r 

axis. The next rotation is about the once transformed X ' = X axis until 

n r 

the coordinates are aligned as follows: 0 about Xj’, 270 degrees about X 2 ', 

180 degrees about X 3 ', and 90 degrees about X 4 '. The results are sum- 
marized as follows in the form of equation ( 10) for each CMG. 


X, = 


C/3 0 


0 1 


S/3 0 


-S/3 

0 

C/3 1 


X 


n = 1 


( 18 ) 


C/3 S/3 0 


X, 


0 0 


X , n = 2 
r 


< 19) 


S/3 -C/3 0 
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C/3 o 


X, = 0 -1 O X , n = 3 


S/3 o -C/3 


C/3 -S/3 0 


0 0 -IX , n = 4 

r 


S/3 C/3 0 


As given by equation (9) , the transformation between spacecraft 
reference and CMG constant momentum coordinates is 


0 0 


X = 0 

c 


Ca Sen X , n = 1, 2, 3, 4 
c c n 


-So Ca 
c c 


Equation (11) is obtained by substituting equations (18) , (19) , (20) , and (21) 
into equation (22) and carrying out the matrix multiplications with c = n . 

For four-skewed CMGs, the transformations between body and CMG 
constant momentum axes are summarized below. 

CMG Number 1 


X t = G. X , G. 

1 lr r lr 


C/3 

0 

-S/3 

S/3Saq 

Caq 

C/3Sa 

S/3Ca 1 

-Saq 

C/3Ca 
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CMG Number 2 




C/3 

S/3 

0 


*2 = G X , 
2r r 

G 2r = 

S/3 So 2 

-C/3Sq 2 

Ca 2 

(24) 



S/3CO, 

-C/3Co o 

- Scio 


CMG Number 3 


C/3 

0 

S p 


X 3 = G„ X , 
15 3r r 

G 3r ' 

S/3So 3 

-Co 3 

i 

-C,3So 3 

(25) 



S/3Co 3 

SO 3 

-C,iCo 3 _ 


CMG Number 4 


" C/3 

-Sp 

0 


X 4 = G X , 
4 4r r 

°4r 

S/3 So 4 

C/3So 4 

-Co 4 

(26) 



_ S/3Co 4 

Cp'Co 4 

So 4 


Utilizing body to CMG transformations, the momentum for each CMG 

can be written in body c 

oordinates [equation (14)] as 


n 

1—i 

II 

l-c 

(S/3So< i 

1 i 

+ Co j j 

i J r 

+ C/3So j 

M \ 


^2 = ^2 j 2 = h 2 

(S/3 So; 2 i 
6 r 

, - C/3So 2 

J r + CO 2 

V f 


^3 ~ h 3 j 3 - h 3 

(S/3 So: 3 i 

o r 

- Co 3 j 
j J r 

- C/3 So 3 

v( 

( 27) 

^4 = h 4 j 4 - h 4 

(S/3So 4 i 
4 r 

, + C/3So 4 j - Co 4 

V 
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The total CMG momentum is the vector sum ol all CMG momentum 
vectors [equation (16)]; thus, 

4 

H = 7) h = h i + h j + h k , (28) 

LJ . cxryrar 
c=l J 

where 

h^ = S/3 (hj So^ + h 2 So; 2 + h 3 Sa 3 + h 4 Sct 4 ) , 
hy = h 4 C&i - h 3 Caj + C/3 (h 4 Sa 4 - h 2 Sct 2 ) , 

and 

h = h 2 C& 2 - h 4 Cct 4 + C/3 (h t Sa, - h 3 Sa 3 ) 
z 


As previously stated, the CMG momentum in reference coordinates will be 
used as the basis for momentum management to prevent CMG saturation and 
to make the CMGs operate about their null positions. The components of 
equation (28) are zero when the gimbal angles are zero. However, there 
are other combinations of gimbal angles which also produce a null momentum 
condition. 

Using equations (23) through (26) , the individual CMG torques 
[equations (8) and (15)] are obtained in reference coordinates as follows: 


hi 

h 2 

h 3 

h 4 


a 1 hj (S/3CQ i - So, j + CPCa i , k^) 
a 2 h 2 (S/3 Co; 2 i^ - CfiC® 2 _ Sa 2 k ) 
a 3 h 3 (S^Ca 3 + Sa 3 - C^Ca 3 k r ) 

a 4 h 4 (S/3C a 4 + C/3Ca 4 j + Sa’ 4 k ) 


(29) 
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The total CMG torque, equation ( 17) , is obtained by summing the contribu- 
tions from each CMG: 


4 



c - 1 


h 

c 


hi + h j + h k 
x r y r z r 


(30) 


where 


h hi Cctj + cVg h .2 C & 2 oi 2 h 3 Ccv 3 t- h^ 

• • • • • 

h y - h 4 Sa t - a 2 h 2 C/3 Ca 2 + a 3 h 3 Scv 3 + a 4 h. C/SCa 4 , 

and 

• • • • . 

h ^ 1 hj C/^Cci j - ct'o h 3 So.' 2 — ci 3 h 3 Ci^Cc*. 3 + cv ^ Set ^ « 


Equation ( 30) can be arranged in the vector matrix form 


- • 


- 




— 

h 

X 


h j Sj3Ci\ j 

h 2 S 1/3 C a 2 

h*> SpC 0 .> 

hj *vtV. j 


“ 1 

h 

y 


-hi SG j 

-h 2 Cp’Cft 2 

h 3 So 3 

h 4 C/3 Ccv 4 


a 2 

1 

tr • 

N 

1 


hj CpCa i 

— h^ Sg 2 

-h 3 CPCa 

h 4 Sc. 4 


a ’ 0 

• 1 






J 

_ a 4 _ 


(31) 


In compacted notation, equation (31) is written as 


H 


C a 


(32) 
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where C is a 3 by 4 matrix denoted as the CMG torque matrix and H and 
a are column vectors. Notice that the columns of C are vectors directed 
along each CMG torque axis, . Since there are four torque vectors, the 

columns are linearly dependent. In the foregoing sections, the momentum 
and torque potentials for the baseline four-skewed CMG configuration have 
been developed relative to the spacecraft reference axes. The next steps 
are to select a skew angle and to examine several candidate steering laws. 


E. Skew Angle and Momentum Capacity 

The foregoing equations have been derived without selecting a specific 
value for the CMG skew angle P , which has been assumed to be equal for all 
CMGs. Several factors enter into the selection of P : (1) momentum capac- 
ity per axis and total momentum envelope, (2) control torque capability 
around the null position, (3) alignment of each gimbal axis to provide the 
independent degrees of freedom required for three-axis control. When one 
CMG has failed, the remaining three CMGs must be able to control the vehicle 
without degrading performance. With this in mind, a skew angle of 45 degrees 
would provide the greatest angular distance between gimbal rate vectors and 
between reference and gimbal axes. The CMG system would, therefore, 
provide the best operational capability with one CMG out. If the skew angle 
were 90 degrees, control torques could be attained about each reference axis 
but the X r axis would have twice the momentum storage capacity as the other 

two axes. Moreover, with one CMG out, severe cross coupling would result 

on the X axis by trying to command only a Y or Z torque. For 
r r r 

example, if fq were out, with P equal 90 degrees, only CMG Number 3 
could produce a axis torque, but that torque could not be produced with- 
out also torquing the X and Z axes. 

r r 

The skew angle could be selected to give equal torque capability per 
axis near the CMG null position. By setting the gimbal rates to some pre- 
determined upper limit (depending on the CMG torque motor characteristics) 
and setting the sign to give maximum torque per axis, equation (31) at the 
null position reduces to the following equations: 


h ( max) 
x 


h ( max) 

y 


h (max) 
z 


4h S/3 a 

2h Cp a 
2h C P a 


( 33 ) 
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Equating maximum torque components produces 


tan (13) = 0. 5 


( 34 ) 


A skew angle of 26.6 degrees, therefore, provides equal torque per axis 
capability near the CMG null position. However, the momentum envelope is 
not symmetric and, as the gimbal angles vary, the torque capability per axis 
does not stay equal. Since the gimbal angles may become rather large if 
momentum is dumped infrequently, equal torque per axis at the CMG null 
does not appear to be a good criterion for selecting the skew angle. 


A more logical approach is to select the skew angle so that the CMG 
momentum envelope is spherical, that is, equal momentum capacity per axis. 
By setting the gimbal angles to values which produce maximum momentum per 
reference axis, equation (28) reduces to the following equations: 


h (max) 

x 

h ( max) 

y 

h (max) 

z 


4h S/3 

2h ( 1 + CP) 
2h (1 + C/3) 


( 35 ) 


Equating maximum momentum components produces 


2 S(3 - 1 + CP 


( 36 ) 


By squaring each side and eliminating S 2 /3 by trigonometric identity, the 
following quadratic equation is obtained: 


5 C 2 /3 + 2 C/3 -3 = 0 


( 37 ) 


The solution of equation (37) gives a skew angle of 180 or 53. 1 degrees. 
However, 180 degrees is a false solution since the axis momentum w'ould 

be zero. Table 1 gives the maximum momentum capability per axis for several 
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TABLE 1. MAXIMUM MOMENTUM CAPABILITY FOR 
FOUR-SKEWED CMGs (ft-lb-sec) 
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skew angles and CMG momentum values. With a skew angle of 53. 1 degrees 

and a unit CMG momentum of 250 ft-lb-sec, each axis has a CMG momentum 

potential of 800 ft-lb-sec for control purposes. A skew angle of 28. 1 degrees 

gives twice as much momentum on the Y and Z axes as on the X axis, 

r r r 

whereas 45 degrees gives 150 ft-lb-sec less on the X axis as compared to 

the transverse axis. Based on the foregoing analysis, a skew angle of either 
45 or 53. 1 degrees is recommended for HEAO-C. 

An analog computer program has been developed 1 to determine the 
maximum momentum surface that a particular set of four single-degree-of- 
freedom, skewed CMGs can generate. The CMGs are mounted so that their 
momentum vectors always lie in the planes of the faces of a pyramid, as 
shown in Figure 3. As each CMG is gimbaled, its momentum vector will 
rotate in the plane. No gimbal position or rate limits are placed on the CMGs, 
and it is assumed that they have equal momentum. The angle of inclination 
of the faces of the pyramid (/?) may be varied from 0 to 90 degrees. 

There is always some total momentum vector H which is the vector 
sum of the individual CMG momentum vectors. If the four CMGs were caused 
to rotate in a random fashion, the locus of the tip of the total momentum vector 


1. W. J. Weiler, PD-DO-ES, MSFC, contributed the material on CMG 
momentum envelopes. 
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would describe a solid. The boundary of this solid is the desired maximum 
momentum surface, or momentum envelope. This envelope is a function of 
the physical system and is independent of the control law used to command 
the CMGs. However, a control law is a necessary part of the scheme used 
to generate plots of the momentum envelope with the maximum contribution 
steering law presently being utilized. The accuracy of this control law 
determines the conformance of the representation to the actual envelope. 


The program commands a total momentum vector of greater magnitude 

than the system can produce. This commanded vector remains fixed in length 

and follows a prescribed pattern in direction. It begins pointing up the 

+ X-axis. It then increments through a fixed angle in the X -Y plane toward 

r r 

the Y^-axis and then revolves about the X^-axis. It continues incrementing 

and rotating until it reaches the -X -axis. The control law causes the indi- 

r 

vidual CMGs to rotate making the total actual momentum vector follow the 
commanded total momentum vector. The rectangular components of the 
actual total momentum vector are plotted by an X-Y plotter to obtain various 
views of the locus of its tip. Ideally, the actual vector 1 would follow the com- 
manded vector exactly in direction, and would maintain the greatest length 
possible in every direction. Actually, due to sensitivity points and singularity 
points of the control law, there is some deviation in parallelism of the actual 
vector to that commanded in some regions, especially when the commanded 
vector becomes nearly parallel with one of the GMG gimbal axes. In this 
case the other CMGs must provide all the momentum in that direction and 
also cancel out the CMG whose momentum is perpendicular to that direction. 


One or more CMGs may be failed by setting its momentum to zero. 

No modification to the control law is required when the CMGs are failed. 
Figures 6 and 7 show profiles of the momentum envelopes for (3 = 53.1 
degrees with all CMGs operational and with one CMG failed. The figures are 
scaled in terms of normalized momentum where one major graph division 
represents one H, the momentum of one CMG. Much distortion of the 
surface is observed when a CMG is failed. The white areas centered about 
the gimbal axis should be interpreted as depressions in the surface, not as 
holes extending through the solid. The absence of contours in the regions 
is mainly due to deviation of the actual vector from the commanded because of 
control law sensitivity points. Depressions do exist there and have been 
verified by digital computer simulations. 
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Figure 6. Four-skeweci CMG maximum momentum envelope. 
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Figure 7. Maximum momentum envelope with CMG number 4 failed. 

When the vehicle is rotated, as in a normal 90 degree maneuver, 
momentum is interchanged between vehicle axes. Simulations indicate that 
environmental torques cause CMG momentum to accumulate in a bias direction. 
These torques depend upon orbital and attitude hold conditions. By the proper 
maneuver, the accumulated momentum can be transferred to any other axis. 
Since HEAO has to maneuver often, the phenomenon of momentum transfer 
by maneuvering makes a near-spherical CMG momentum profile highly desir- 
able. A skew angle /3 of 53.1 degrees will produce die near-spherical 
momentum envelope shown in Figure 6; therefore, it has been recommended 
for HEAO. There are slight indentures at each CMG gimbal axis because a 
CMG cannot contribute any momentum along its gimbal axes. At the inden- 
tures, I Hi is about 675 ft-lb-sec for 4-250 ft-lb-sec CMGs. With one CMG 
failed, [H| at the indenture is about 368 ft-lb-sec for diree 250 ft-lb-sec 
CMGs. It should be noted that doubling die j H I per CMG would also double 
the size of the momentum envelope. 
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SECTION III. CMG STEERING LAWS 
A. Introduction 

When the CMG gimbal angles are moved, a corresponding change 
occurs in the momentum relative to the body axes. By definition, torque is 
the time rate of change of angular momentum. Therefore, the spacecraft is 
acted upon by a torque when the CMG gimbal angles are changing. In general, 
a spacecraft control law is derived as a linear combination of sensor outputs 
such as rate gyros, sun sensors, star trackers, etc., which have been 
weighted by a constant gain on each output. The gains are selected to give 
the desired vehicle response and stability characteristics. Ideal control 
would be obtained if the torque called for by the vehicle control law could be 
produced by the CMGs. The control law is typically derived with respect to 
the spacecraft reference axes and may be written in the following general 
vector form: 


T — T i+T j+T k 
c cx r cy r cz r 


(38) 


The basic objective is the derivation of a CMG gimbal control law providing 
some approximation of the torque specified by the spacecraft control law. 

The standard approach is to equate the total change in CMG angular 

momentum H , equation ( 32) , to the desired torque T . The vector com- 

c 

ponents are equated with a negative sign and arranged in the following vector 
matrix form: 


~T 

cx 




• 





a 2 

T 

cy 

= 

3 x n 
Matrix 


* 

T 

cz 



i 

• 

1 a 

n 


( 39 ) 
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The 3 by n matrix must be inverted to obtain a general solution for the CMG 
gimbal rate commands. Several problems are obvious: ( 1) With lour or 
more CMGs, the linear system is underdetermined, meaning that when the 
equations are consistent there is an infinite number of solutions (there are 
only three equations but n unknowns) : (2) For some gimbal angles, the 
system is known to be inconsistent and not all sets of gimbal angles producing 
inconsistency have been determined ( lor some gimbal angle combinations no 
solution exists) ; (3) These considerations and the algebra involved make a 
general solution almost impossible without resorting to a digital computer. 

A second approach is to restrict the range of the CMG gimbals and 
assume small deviations from the CMG null positions, in such a case, small 
angle approximations are used, sin n - <v and cos o 1 , and the equations 
are linearized. The gimbal rates are solved so that cross coupling between 
axes is eliminated. The resultant solution, however, is valid only for small 
CMG angular excursions from the null positions. At this point in the CMu 
control system design, each individual designer will have or devise his own 
method for selecting a CMG steering law. Several candidate steering laws 
were derived and compared on the basis of their effectiveness in producing 
the desired actuatoi’ response as well as their complexity in implementing 
each scheme. Each steering law was derived for the four- skewed CMG mo- 
mentum exchange system which has been baselined for HEAO-C. 

For the four-skewed CMG configuration there are three equations 
(components of the CMG torque vector) and tour unknowns 1 tour CMG gimbal 
rates) . To obtain an exact solution, a constraint equation or relation between 
the unknowns is needed. For each constraint or assumption that is made, a 
solution will be obtained for the gimbal rates. Whatever the solution, it is 
referred to as the CMG steering law. The following are several candidate 
steering laws that were considered: 

1. Constant gain. 

2. Maximum contribution. 

3. Psuedo inverse. 

4. The Bendix summation ol three-gimbal inverses. 

5. The BECO momentum vector distribution. 

6. The G. E. transpose with torque feedback. 



Because of the large quantity of data, simulation results for each steering law 
are included in Appendix C , including magnetic momentum management for 
the CMGs. 


B. A Constant Gain 

A constant gain steering law can be derived by assuming that each CMG 
will operate about its null position and that torque must be generated about any 
vehicle axis. The CMG torque in body axes is equated to the desired control 
torque to obtain the following variation of equation (39): 


- T cx - hS/i(cq Ca^+ & 2 Ca 2 + &3 Ca 3 + <a 4 Ccv 4 ) 

= h( -Oj Soj - Oi 2 C /3 C a 2 + ® 3 So , 3 + 0^4 C(i Ca 4 ) 
" T CZ = h( " 1 C(i Cftl - Sa 2 - a 3 c /3 Cq 3 + a 4 Sa 4 ) 


(40) 


Assuming small gimbal angles, Ca.= 1 and Sq = 0 , equations (40) 
reduce to 


- T Cx= hS ^<«l+ «2+ 0,+ Q 4 ) 

-T = hC /3 ( a 4 _ o ,) 
cy * l 

-T =hC/3(o,-aJ 
cz 1 6 


(41) 


For either torque or momentum capability, CMGs number 2 and number 4 
dominate the Y-axis, and CMGs number 1 and number 3 dominate the Z-axis. 
Any CMG can be used to produce torque on the X-axis. Since there are four 
unknowns in the gimbal rates but only three equations, the solution for gimbal 
rates as a function of desired control torques is not unique. Some criterion 
or constraint between the gimbal angles is needed for a unique solution. 
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Assume that only X-axis torque is wanted and that it must be produced without 
introducing torque on the transverse axis. By setting a 3 = Oj and a 4 = d 2 
both the Y and Z torque components are zero and, by setting cv 2 ^ aq 
the X-axis torque attains a maximum value of 


-T - 4 hS pet, 
cx 


(42) 


Solving equation ( 42) for the gimbal rates produces 


a 


xi 


T 

cx 
4 hS P 


i = 1, 2, 3, 4 


< 43) 


where the subscript x represents the X-axis solution. 

Similarity, the Y-axis torque is maximized and the transverse torques 
are zeroed by setting a 2 = -oq and oq = a-, = 0 . The Y-axis torque com- 
ponent is 


-T = hC/i( 2 oq) 
cy 

The corresponding gimbal rate solutions are 


( 44) 


a 


yi 


= « , = 0 


o 

y^ 



= T / 2 hC 13 
cy 


(45) 


The Z-axis torque is maximized by setting Q :j - - a 1 and a ■> = a 4 = 0 . 
Using these values, both the X and Y torque components are zero and 


-T - hC £ ( 2 a , ) 

n y 1 


( 4C) 
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The gimbal rate solutions are 


a = -a = -T /2 hC/3 
zl z3 cz 


“z2““ Z 4 =0 


( 47 ) 


The constant gain steering law is obtained by summing up equations (43) 
(45) , and (47) as follows: 


a=a +a +a = -T /4 hS/3 - T /2 hCj3 
1 xl yl zl cx cz 


a 0 = a 0 + a + a . , = -T /4 hS/3 + T /2 hC/3 
2 x2 y2 z2 cx cy 


a = a + a + a , = -T /4 hSi3 + T /2 hC P 
3 x3 y3 z3 cx cz 


oi - a .+ a + a = -T /4 hS/3 - T /2 hC/3 
4 x4 y4 z4 cx cy 


By defining constants 


(48) 


K A = 


k b = 


4 hS/3 


1 


2 hC/3 


( 49 ) 


equations (48) can be written as 



«3 = K. T - K T 

A cx B cz 


«.i K T + K T 

A cx B cy 


( 50) 
( Cont'd) 


Equations ( 50) iilustrato the constant gain CMG steering law that was first 
used in the HEAO-C CMG performance simulation studies. At a later date, 
the maximum contribution steering law was derived, which reduces to the 
same constant gain steering law by assuming small gimbal angles and 


linearizini 


*&• For a skew angle of 55. 1 degrees and 250 It-lb-see C.MGs. the 
gain constants are 


K, = -0.00125047 
A 

K = -0. 00533111 
B 


< 51) 


The constant gain steering law is most simple to implement and could 
easily be simulated on an analog computer. However, it is valid only tor 
small gimbal angles. In the digital simulations, excellent performance was 
obtained if the gimbal angles were less than ±4 5 degrees. For many cases, 
satisfactory performance was obtained with gimbal angles up to ±80 degrees, 
Conti ol was lost it the angles exceeded ±90 degrees. With continuous momen- 
tum, dump using magnetic coils, for example, the gimbal angles stay small and 
the constant gain steering law meets all IIHAO-C requirements. However, 
with periodic momentum dump using RCS thrusters, less than half the avail- 
able momentum can be used belore the gimbal angles exceed their linear 
operating range. For the baseline HFAO-C configuration with four 250 l’t- lb- 
sec CMGs. momentum would have to be dumped each one- hall orbit under 
worst-case environmental torque conditions. As an alternative, a more 
general type constant gain steering law with periodic gain switching could be 
defined to permit better utilization of the total momentum capacity. The latter 
approach was not pursued during this study. 


With one CMG tailed, the CMG null position must be redefined and a 
constant gain steering law derived which is valid about the new null position. 
The dump frequency would have to be increased to about four times per orbit 
assuming worst-case environmental effects. Hence, for each CMG failure, 
a contingency steering law' must be oofined witn a new’ CMG null position. 
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With continuous momentum dump with electromagnets, any CMG can be failed 
and a new null automatically found without reprogramming. However, vehicle 
maneuverability is restricted by the linear operating range of the steering 
law. So even with continuous momentum dunip, the constant gain steering 
law should be redefined for each failure mode. 

Assume that CMG number 3 has failed. With the failed CMG eliminated, 
the momentum becomes 


h = hSP(Sa i + So 2 + So^) \ 

X 1 

h y = MCffj - C/3 Sq 2 + CpSai) \ . (52) 

h^ = h(C/lSo'i + Cct 2 - C 04 ) J 

One new null position can be found by setting 0 ^ = 0 and solving for a 2 and 
a 4 . With q ?4 = 0 the X-axis equation gives 02 = -0-4 , which also produces 
zero in the Z-axis. The Y-axis momentum set to zero gives 


0= h(l - CfiSa 2 + C/lSo 4 ) 

With a 2 - -a 4 , equation (56) reduces to 


With (3 = 53. 1 degrees, the solution for a new null position is 
a j = 0 deg 
0! 2 = 56. 4 deg 
«4 = -56. 4 deg 



(53) 


(54) 


( 55 ) 
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The CMG torque potential with CMG number 3 out is [equation (40) with 
So?! = 0, Coq = 1, Ca 2 - Cft 4 = 0. 5534, Sa 2 = 0. 8329, and So 4 = -0.8329] 


ii = -T = hS/K ft i + 0.5534 a, + 0.5534 ft,) 
x cx 


h - -T = hCj8( 0.5534) (a , -ft,) 

y cy 


h = -T = h{Cfta , - 0. 8329 a, - 0. 8329 ft .) 
z cz 


( 50) 


Since there are three equations and three unknowns, an exact solution can be 
found, assuming the equations are consistent. In vector matrix form, equa- 
tions (56) become 


1 

0. 5534 

0. 5534 

( « i ) 

0 

-1 

1 

1 a 2 1 

Cft 

-0. 8329 

-0. 8329 

(ft 4 ) 


-T /hSft 
cx 

-T / 0. 5534 hC ft 
cy 


-T /h 

cz 


< 57) 


The determinant of the matrix, A , for ft = 53.1 degrees is 


A = 2(0. 8329 • 0. 5534 Cft) - 2. 3300 


< 58) 


Since the determinant is not zero, a solution is found by using Cramer’s rule. 
With h - 250 ft-lb-sec, the constant gain steering law is given by the following 
equations: 


ft, = -0. 003570 T 

cx 


ft, = -0. 001289 T 

cx 


ft.- -0. 001289 T 

cx 


0. 0019 T 

cz 


0. 00002 T - 

cy 


0. 00602 T + 
cy 


0. 001716 
0. 001716 


T 

cz 


T 

cz 


I 59) 
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Substituting the steering law, equation ( 59) , into the CMG torque equations, 
equations ( 56) , the CMG torque per unit command torque is 


h /T 

x 


cx 


h /T 

y 

h /T 

z 


cy 


cz 



(60) 


Thus, with h 3 = 0 , the gain through the CMG system on any axis is unity, 
and the steering law solution is verified. A similar procedure would be 
required to obtain a solution for any other CMG out. 


C. The MSFC Maximum Contribution 

A second and somewhat novel approach is to command each CMG 
separately based on its ability to contribute to the desired control torque. 

The criteria are to consider each CMG independently and to command its 
gimbal rate so that as much as possible of the desired control torque is 
produced. If no part of the control torque can be produced, the gimbal posi- 
tion is not moved. Since each CMG can produce a torque only about its torque 

axis as defined by k , the desired control' torque T will be projected into 

c c 

the CMG constant momentum coordinates and used to command the gimbal rate. 
Using the transformation G , the desired torque, as defined by the control 

law, in CMG coordinates is as follows: 


T = t i + 1 j +t k 
c cx c cy c cz c 

where 


t 

cx 

c 

= §11 

T 

cx 

c 

+ 8 12 

T 

cy 

c 

+ §13 

T 

cz 

> 

t 

c 

§21 

T 

c 

^ 822 

T 

+ §23° 

T 

J 

cy 

cx 

cy 

cz 



(61) 
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and 


t 

cz 


c 

£31 


T 

cx 



T 

cy 


e 


O’ CO 
^ oo 


T 

cz 


Previously , the torque produced by the cth CMG was defined as 

• » 

h a h k 

c c e c 


By equating elements of 6, and -T, , the CMG gimbal rate command Is 
obtained as 


a 


c 


-t h 

cz c 


-(o'.,.. 


T 

cx 


f>32 


T 

cy 



T 1 h 
cz c 


( 63) 


By defining the angle between T and k as -y 

c c 

is 


the vector scalar product 


T . k 

e c 



T 

cx 


+ 


T 

cy 


+ 




T (>. 

c 


( 64 ) 


therefore, 


a 


-T CV /h 
c c 


( 65) 


Note that equation (64) is the projection of the desired control torque on the 
k c axis, i lat portion of y which is perpendicular to too CMG torque axis 


is given by 



( 66 ) 
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and cannot be obtained by gimbaling the cth CMG at any time. Since T Cy 

c r 

contains basically magnitude information without polarity, the expanded scalar 
product form should be used for the CMG gimbal commands. Each CMG is 
commanded individually, regardless of its angular position, based on its 
ability to contribute to the desired torque vector. For a specified CMG system 
configuration, however, the transformation between body and CMG coordinates 
must be derived and evaluated at each time step. 

If the components of G from equations (23) through (26) are sub- 
stituted into equation (63), with c = 1, 2, 3, 4, the following equations are 
obtained for the gimbal rate commands for four-skewed CMGs: 


“l ~ 


01 2 = 


- T cx - Sa, T cy + CPCu, T_ ) /h, 


cz 


- (seca 2 T cx - C/3Ca 2 T_ - Sa 2 T J h 


cy 


CZ 


« 3 = - (SpCa 3 T cx + Sa-3 - Q3Ca 3 TJ ,/h 


cy 


cz 


(67) 


a 4 = 


<SSCa * T cx + CKa > T cy + Sc > ' h 4 


CZ 


Assuming equal momentum per CMG, the steering law can be arranged in the 
following vector matrix form: 


a = A T /h 
c 


( 68 ) 


As previously shown by equation ( 32) , the CMG torque with equal momentum 
per CMG is 


H = h C o' 


(69) 


By comparing the elements ot A and C , the matrix A equals the negative 
transpose of C ; that is, 


A - -C* 


( 70 ) 
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Hence, the CMG steering law shown in equation (67) is equivalent to a trans- 
pose type steering law. The important characteristics of this type steering 
law are no mathematical (computation or algorithmic) singularities, no 
matrices to invert, and easy implementation of the equations for the CMG 
gimbal rates. Its undesirable characteristic is that the control system gain 
through the CMGs is not constant. Control authority about each vehicle axis 
varies as the gimbal angles are moved from their null position since at each 
instant of time, each CMG is maximizing its contribution to the desired 
control torque. 

Assuming small gimbal angles, the CMG torque per vehicle axis, 
equation (30), reduces to the following equations: 

• • • • • 

h h S/j ( O' i + ot 2 + Cl ^ -t- 0 : 4 ) 

» • o 

h = h C/3 (-or, + G’d 

y 

• ♦ • 

h z h C/3 ( O! j — G g) 

Substitution of equations (67) into equations (7l) results in the following- 
equations : 

h -4 S'V T 

X cx 

h -2 C 2 ii T 

y cy 


h -2 C 2 /i T 
z cz 

To obtain unit gain through the CMG system with the CMGs at their null 
position, either the desired torque components, T^ , or the steering law- 

must be normalized by dividing by the appropriate sine and cosine function of 

the skew angle. If one chooses to normalize the steering law, those terms 

with T are divided by 4S 2 (i and those terms with T or T are 
cx cy cz 

divided by 2C 2 jS . The maximum contribution steering law ensues from the 

normalization of equations (67) as follows: 





Of = 

(-Ca 1 /4h 

S/3) 

T + 

(S <V2h C 2 /3) 

T 

• 



cx 

cy 

a 2 = 

(-Ca 2 /4h 

S/3) 

T + 

(Ca 2 /2h C/3) 

T 




cx 

cy 

a 3 = 

(-Ca 3 /4h 

S/3) 

T - 

(Sa 3 /2h C 2 /3) 

T 




cx 

cy 

“4 = 

(-Ca 4 /4h 

S/3) 

T - 

(Ca 4 /2h C/3) 

T 




cx 

cy 


cz 


cz 


• 3 / 2 h c/5) T 


( 73 ) 


cz 


cz 


Once specified, the skew angle is constant. By making the following defini- 
tions. 


K & = -l/4h S/3 


\ = -l/2h C/3 

V C/? 


(74) 


K = 
C 


the maximum contribution steering law can be written as follows: 

K Ca,T - K Sccj T +K Ca, T 

a cx c cy b 1 c 

K C a 2 T - K C a 2 T - K Sa, T 
a cx b cy c 1 c 


Cti 


an = 


cz 


cz 


a, = 


OtA = 


a « cx c 0 'ey ' ^ C “ 3 T C z 


K_ Ca 3 T_ + Sa 3 T 

K Ca 4 T + K C« 4 T + K Sa, T 
a cx b cy c 4 i 


(75) 


cz 


The maximum contribution CMG steering law has been used extensively in the 
HEAO-C simulations with four-skewed CMGs. Very satisfactory pointing 
performance was obtained in the cases simulated. Although there are no 
computational singularities in the maximum contribution steering law, there 
are gimbal positions for which a component of the desired torque cannot be 
produced. This condition for the maximum contribution law is ghu. For 
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example, assuming CMG saturation in the X-axis direction with all gimbal 
angles at 90 degrees, no torque can be produced in the X-axis direction, 
h = 0 . Concurrently, the desired torque component T cannot drive the 

X C X 

CMG gimbals through the steering law; the coefficients of T in equation 

cx 

(75) are zero. In general, however, ghu conditions are dependent on the 

CMG system configuration and are not necessarily associated with the ability 

to command through the steering law. Consider the ghu condition cv, = 90 

degrees, as =.0 degrees, cv 3 = -90 degrees, and n 4 = 0 degrees. The Z-axis 

CMG torque, h , is zero but the CMGs are not saturated. Also the Z-axis 
z 

command T cannot be fed through the steering law. The implication is 

that for the transpose type steering law, ghu conditions also correspond to 
gimbal positions which prevent a commanded torque component from driving 
the CMG gimbal angle rates. 


D. The Pseudo Inverse 

By equating the desired control torque, T^ , [equation (38)1 to the 

CMG torque, H, [ equation (30 )] the vector-matrix equation, equation (32) , 
becomes 



C (hr* ) 


(76) 


where C has been normalized by factoring out h = h , i = 1, 2, 3, 4* The 

matrix C is a 3 by 4 whose inverse must be obtained to solve the gimbal 
angular rates as functions of the desired control torque. However, since 
there are four unknowns but only three equations, a general solution, if one 
exists, is not unique. One, therefore, must resort to a pseudo inverse. 

Reference 8 gives the general conditions and theory for finding a 
pseudo inverse solution, one form of which can be given by 


F = C : (C C*) 1 


(77) 
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provided that the matrix C is of rank three or, what is the same, that the 
instantaneous torque vectors are not coplanar. The steering law then becomes 

(ha) = -C* (C C*)" 1 T = -F T , (78) 

C c 

where superscript asterisk represents the transpose of a matrix and minus 
one represents the general inverse of a matrix. For the baseline CMG con- 
figuration, the elements of C have been defined by equation (31 ) which must 
also be normalized by factoring out h = h , i = 1, 2, 3, 4 . To obtain the 

inverse, the determinant of CC* must be calculated. Those gimbal angle 
combinations which make the determinant go to zero and cause program 
divergence are denoted as singularities. For the pseudo inverse, singularities 
are synonymous with ghu conditions. However, other steering laws can have 
singularities without a corresponding ghu. Since mathematical manipulations 
required to obtain the pseudo inverse are too complicated to perform without 
a digital computer, the elements of equation (77) will be developed only to 
the extent required for calculations. Most digital subroutines for matrix 
inversion accept the elements of the matrix to be inverted and give as outputs 
the elements of the inverse matrix. Let 


D = C C* 


(79) 


and 


E = D 1 


(80) 


The normalized elements of C are obtained from equation (31) from which 
the elements of D are calculated as shown in the following: 
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Du = 

C ll 2 + C 12 2 + C 13 2 + C 14 2 


Dl 2 - 

Cn C 21 + c 12 c 22 + C 13 C 23 f C !.j 

W, 

D 13 : “ 

C 11 C 31 + C 12 C 32 ! C 13 c 33 + C 14 

C 34 

D 21 ' 

D 12 


D 22 : 

C 2 1 2 + c 22 2 -v c 23 2 c 24 2 


D 2 3 ; 

C 21 C 31 + c 22 C 32 + c 23 c 33 + c , . 

c :;i 

D 31 

Dl 3 


d 32 -- 

D 23 


d 33 

C’J-r Cj-r C 33 2 -r c 54 2 . 
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As indicated above, the matrix D is skew symmetric. 


The elements oi D are inputs to a digital matrix inversion routine 
to obtain the elements oi E as outputs. Multiplying the elements oi E by 
C gives the elements of F , a 4 by 3 matrix, as shown below. 


F 11 C 11 E 11 + C 2 i F 21 + C 31 


F12 

Cll 

E 12 

+ c 21 

E22 + C31 
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F 13 
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C 13 
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II 
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+ C 23 

E 22 - C33 
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(82) 

( cont'd) 


F 33 - 

C 13 E 13 + C 23 

E 23 + E 33 E 33 

II 

tH 

C 14 E 11 + C 24 

E 21 + C 34 E 31 

ii 

CM 

Cl4 E 12 + C 24 

E 22 + E 34 E 32 

OS 

II 

c t4 e 13 + c 24 

E 23 + ^34 E 33 


Utilizing the elements of the pseudo inverse matrix, the CMG steering law is 
obtained by expanding equation (78) as follows: 


a< = 


a, = 


a, = 


a a = 


= -(Fn 
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E 12 
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cx 


cy 

= ~< E 21 
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E 22 
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cy 

= -( E 31 

T 
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23 i cz // u 2 


' 33 


CZ 


CZ 


(83) 


As noted previously, the columns of C are CMG torque vectors for each 
CMG. It can be rigorously proven that when any three are colinear the 
determinant of C C* goes to zero, producing a singularity in the steering 
law. There is a large number of gimbal angle combinations that can produce 
singular conditions. However, in digital simulations, the only singular 
conditions which prevented proper operation of the steering law were those 
which also corresponded to CMG saturation. 


In most cases when an internal singularity was approached, a small 
pointing error was observed. However, the system would recover and oper- 
ate satisfactory until CMG saturation was reached. The cyclic nature of 
environmental torques prevented ghu at the singularities. But when constant 
torques were commanded, internal singularities could always be encountered 
with subsequent loss of control. Only about 50 to 60 percent of the momentum 
envelope is usable without any possible singularities. More research is 
needed to fully understand the singularity and ghu conditions associated with 
single gimbal CMGs and to develop possible avoidance schemes. 
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When a CMG fails, it is acceptable to set the failed gyro elements to 
zero. The failed CMG must be identified and the column corresponding to the 
failed CMG set to zero. The pseudo inverse routine need not be repro- 
grammed. With one CMG out, however, C reduces to a 3 by 3 matrix whose 
inverse can be obtained without resorting to the pseudo inverse procedure. 

The advantage of the pseudo inverse steering law is that most of the time the 
exact torque needed for attitude error correction can be obtained through the 
CMGs without any cross coupling. Possible disadvantages are the complexity 
of implementing the pseudo inverse matrix inversion routine and the require- 
ment for detecting and compensating for CMG failures. An onboard digital 
computer would be required to implement the pseudo inverse steering law. 


E. The Bendix Three-Gimbal Inverse 

In the foregoing sections, the basic characteristic of the four-skewed 
CMG system is that there are more control variables, gimbal rates, than 
there are basic relations, torque equations, between the variables. As 
previously shown, the three components of the CMG torque vector provides 
three equations which can be arranged in the vector matrix form 


H - C (h o') 


(*n 


where C is a 3 bv 4 torque matrix whose columns correspond to unit vectors 
directed along each individual CMG torque axis. From the basic definition 
of the CMG reference coordinate systems defined in equation (8), the indi- 
vidual CMG torque is 


h = h a k , c — 1 , 2, 3, 4 
c c c c 


(85) 


where k is a unit vector in the CMG coordinate system. The vector compo- 
nents of k in the vehicle reference coordinate system form the elements of 
c 

the cth column of the C matrix. In vehicle coordinates, let k , c - 1, 2, 3, 4, 

be the column vector corresponding to the torque vector of the cth CMG, then 
equation ( 84 ) can be reqritten as 
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H = [k, k 2 k 3 k 4 ] (hS) 


( 86 ) 


After replacing the CMG torque, H , with the desired control torque compo- 
nents, T c , the objective is to solve the torque equations for the CMG gimbal 

rates such that the CMG system generates the exact control torque. However, 
there are an infinite number of gimbal rate combinations which will satisfy 
the torque equations. By specifying some subsidiary condition between control 
variables, an equation between the variables can be obtained which can be 
utilized to eliminate one oi the variables from the torque equations, hence 
reducing the torque matrix to a 3 by 3 matrix which will have a unique inverse 
if the columns are linearity independent. The constraint equation must 
necessarily be based on some preconception of what comprises a desirable 
CMG system state or response characteristic. 


In the absence of a universally accepted subsidiary condition between 
the CMGs, Bendix has proposed "the three-gimbal inverse" steering law 1 6] 
for use on HEAO. If there were only three CMGs, the torque matrix would 
reduce to a 3 by 3 matrix by deleting the column corresponding to the deleted 
CMG. In this case, a unique solution exists for the three gimbal angles, 
assuming that the determinant of the 3 by 3 matrix is not zero. In the Bendix 
scheme, the CMGs are grouped into sets of three and the desired control 
torque is apportioned to each set. Each set of three CMGs is required to 
deliver its apportioned part of the desired control torque. Then, the corre- 
sponding CMG gimbal rate commands are obtained by inverting each 3 by 3 


matrix and summing the results from each set. For the four-skewed CMG 
configuration there are four possible sets of three CMGs which result in the 
following equations: 

H — l k 2 kg kj 

th.ij , 

CO 

of 

II 

j 

3? 

CO 

rX 

II 

•IX 

lh i “i> • 

i = 1, 3, 4 | 


H = [ kj k 2 kj] 

(1 W • 

i = 1, 2, 4 | 

^ • (87) 

H = [ kj k 2 k 3 J 

l VV ■ 

i = 1, 2, 3 


In each equation, let 

A c ' C= *■ 

2, 3, 4, be the torque matrix corresponding 


to the set ot three CMGs with the cth CMG deleted lrom the 3 by 4 torque 
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matrix C. The inverse ot each A exists ii the determinant of A is not 

i c c 

zero. Let A denote the inverse of A and a denote the eorrespondimr 
c c c 

solution for the three gimbal rates based on the cth set of three CAIGs. The 
gimbal rate solutions for the four CMG sets are 


a t - (A,- 1 H)/h 
a L > ( A 2 -1 H) /h 
a 3 ( A ;i -1 H) /h 
a 4 - (A4- 1 H)/h 


CMG no. 1 out 
CMG no. 2 out 
CMG no. 3 out 
CMG no. 4 out 
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Although it is not necessary to prorate the desired torque equally among the 
lour CMG sets, there is no basis for doing otherwise. If the CMG torque 

H is replaced by the desired torque vector T , then conceivably each set 

could deliver the total required torque. To prevent overtorquing. the desired 
torque components are divided by four, that is, apportioned equally between 
the lour sets. With T /4 being substituted into equation ( 85 ) with a negative 

sign, the lour solution sets are obtained and the results for each CMG gimbal 
rate added together to obtain the Bendix three-inverse steering law. A flow 
diagram ol the steering law is shown in Figure 8. The gimbal rate command 
to each CMG is composed ot solutions trom three ot the lour solution sets. 

In the event that one CMG fails, only that set which does not contain the 
failed CMG would be used to obtain the gimbal rate commands. For example, 
it CMG number 3 tails, then all the desired torque would be allotted to the 

A3- solution. In this case. a 3 would give the exact solution needed to 
generate the required control, if and only if k,, k 2 , k 4 are not coplanar. 

When three unit torque vectors are coplanar. the vector box product between 
them is zero. Moreover, the box product is identical to the value of the 
determinant formed by the vectors. When the determinant is zero, no solu- 
tion exists, and the matrix is singular. With CMG number 3 out. the gain 
lactor l/4h would be changed to 1/h and the loops broken which lead to 
^1 * ^2 > and A 4 1 . The resulting solution from A 3 ~* should be, in 

this case, the exact solution. 




To CMG #1 


To CMG #2 


To CMG # 3 


To CMG #4 
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< 89) 

< eont'd) 



where the subscript on the gimbal rate solution sets are used tor identifica- 
tion ol the particular solution. The CMC gimbal rate commands can be 
written in vector matrix form by augmenting each solution with a null (zero) 
row which corresponds to the deleted torque vector: 



That is. each gimbal rate command is composed of the three solution sets in 
which its torque vector appears: 


a i o 12 ‘ 

• ♦ 

e 2 (■ 21 + 

O :j O :j | ~ 

’■ 1 *■ 11 1 


(■ is * e i: 
r. j:i • O ;M 
'■ 32 " e 34 
(■ 42 * *> 43 


( 91) 


The matrix inverses can be incorporated into equation * 90) by adding the null 
row to each inverse matrix to form B a 4 by .9 matrix with its clh row 

zero, for example. 
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0 


0 


0 


a ll 

a 12 

a 13 

a 21 

a 22 

a 23 

a 31 

a 32 

a 33 


( 92 ) 


Utilizing this somewhat superficial matrix definition, the three inverse steer 
ing laws can be written as 



( -1/ 4h) (Bj+ B 2 + B s + B 4 ) 



(93) 


1X31 a ij ’ a i/ ’ a jj 3 » and a jj 4 represent the elements 

A i . A 2 _1 > A 3 _1 i and A 4 1 respectively. Then, from 
(93) , the gimbal commands are 


of the inverse matrices 
equations (92) and 
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a -4 - (-l/4h) 


( a 3i* + a 3i“ + a :ii 3 ) T cx + a 3 | f a 3 |) T 


cy 


+ 1 a.ij 1 + a 3s" 


?4 


(94) 

( cont'd) 


II the inverse exists for all four girnbal sets, the desired torque will 
be obtained. However, if one or more of the A matrices are singular. 

special strategies must be devised to obtain the desired torque. The solution 
set whose determinant is zero could be disregarded and the desired torque 
apportioned to the remaining three sets. Bendix proposes a CMC singularity 
detection and avoidance scheme in their HEAO-A Phase B Final Study Report 

I 6j . I he box product between the column vectors of each A matrix is 

c 

continuously calculated and, when any set value becomes less than some 
specified small value, a biased rate command is applied to one of the three 
CMGs in that set. The remaining CMGs must counteract the torque produced 
by the biased rate, hopefully driving the CMGs away from the singular 
condition. 


Singularity detection is accomplished by continuously monitoring the 
triple scalar product between the column vectors of each torque matrix. A 

The value ol the determinant ol A is identical to the triple scalar product. 

When the determinant of a torque matrix is zero, the three torque vectors 
l represented bv . c 1. 2. :J. 4. in equation ( 87) J are coplanar and that 

particular matrix has no inverse at that instant of lime. Singularities are 
delected by monitoring 
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When the absolute value of any determinant is less than a small positive 
constant, P g , a near singularity has been detected for the cth torque 

matrix. That is, if 


i i a ; ; < p 

c s 


( 96 ) 


the singularity avoidance scheme is invoked. 

comm Jd he f gUlar ; ty ““ S ° heme °° nSiStS °' aPP ‘ ying a blased rate 
command, A a _ c , to any of the three CMG gimbals represented in the matrix 

whose determinant is less than P s . The sign of the bias is opposite the 

r lar ‘‘ y ? f th ® g ‘ mbal rate Just Previous to invoking singularity avoidance. 
Letting „ s be the magnitude of the bias, the bias rate command to the ith 

CMG can be written as 


Acy ir = " Si ^' n ^ ■ tt-l)J 

s ic 


(97) 


" b T 6 lndlcates th ® rale measurement from the previous computational 

^tt' n TT* lhrea CMGs . denoted by j. k, 1 are biased to counter- 
act the bias applied to the ith CMG by defining 


An 


jc 


At* 


kc 


Ao 


lc 


-o- A. 1 k 
s 1 i 


(98) 


where k. is the column vector corresponding to the biased CMG torque 

vector, equation (86) , and is the inverse of the matrix that does not 

contain the ith torque matrix. For example, if iAj, < p > then either 

gimbal 2 3, or 4 may be selected for a bias rate command. Selecting the 
second CMG gimbal to apply a bias rate command yields 
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and 


A <' n - ~ d sign l o 21 ( t-l> J 


< <>;)) 



' J OO; 


Then the bias commands are added to the solution set rate command 
Uon Ji ) . to provide singularity avoidance. For the example given 
gimbal rate commands are 


s. equa- 
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Once initiated, the singularity avoidance- slewing will continue until eitner 
another determinant drops below or the lirst exceeds P . 11 another 

set drops below P g . the bias commands are based on the second singularity 
condition. 11 all determinants exceed 1^ . the- avoidance is discontinued and 
norm aJ opo ration is resumed. 


! n . e P iecedin 8 paragraphs, me liendix three-gnuuai inverse steering 

-aw with singularity detection and avoidance has been derived. 11 each solu * 
o>n set could contribute Us allotted portion ol the desired commanded torque 

and the°dr V, <,bUllno(L I)c '> K ' l,flj ^ »P‘>n the CMC, gimbal positions 

and the desired torque at any instant of time, however, the solution sets 

cannot pi wide equal control authority in the required direction. For example 

-hen any two CMC torque vectors, k . are eolinear. at least turn o, the ' 
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solution set matrices are singular [ equations (87)] and cannot provide their 
apportioned share of the desired control torque. Even if the other two solu- 
tion sets could provide perfect commands, only half the needed torque would 
be produced since each set has been apportioned only one-fourth of the 
required torque commands. Moreover, when the singularity avoidance 
scheme starts working, the bias commands will produce undesirable torque 
components which will interfere with the normal torque needed for control, 
temporarily interrupting the vehicle's pointing accuracy. To help alleviate 
this deficiency, each solution set should be required to produce all the 
desired torque and the individual gimbal commands adjusted by dividing by 
one-fourth. 

In comparing the Bendix three-gimbal inverse steering law with those 
derived in the previous sections, the Bendix law is far more complex then 
any of the other laws. Moreover, the Bendix scheme introduces mathematical 
singularities that are not otherwise present and do not correspond with ghu. 

For example, if any two CMG torque vectors are colinear, two of the solution 
set matrices become singular even if the total CMG system can provide per- 
fect control through the remaining two solution sets. 

Singularity detection and avoidance schemes are required to prevent 
program divergence. When the Bendix singularity avoidance scheme is 
invoked, unwanted torques are introduced which tend to disturb the vehicle's 
pointing performance. Four 3 by 3 matrices must be inverted to obtain the 
Bendix steering law, but only one 3 by 3 inversion is required in the pseudo 
inverse steering law. CMG failures must be identified and the correct solu- 
tion sets deleted for proper failure mode operation with the Bendix steering 
law. Because of the complexity of the Bendix three-inverse steering law 
without corresponding increases in either reliability or performance, it is 
not recommended for use on HE AO. 

Simulations indicate that the Bendix three-gimbal inverse steering law 
produces acceptable vehicle pointing performance. But, without the singularity 
detection and avoidance scheme, only about one- fourth of the available CMG 
momentum could be utilized before encountering a singular condition, after 
which the system diverged and exceeded the required pointing specifications. 

To utilize the Bendix steering law and the total CMG momentum envelope, it 
is absolutely necessary to also use their singularity detection and avoidance 
scheme, which will degrade performance. About four times more computer 
time was required for the Bendix steering law without singularity avoidance 
than for the pseudo inverse steering law. 
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F. The G.E. Transpose with Torque Feedback 

As previously shown by equation (32 ) , the CMG torque is 


H - hC a 


(102) 


where h is the momentum per CMG, C is the normalized 3 by 4 gimbal 

torque matrix, 5 is a 4 by 1 column matrix representing the CMG gimbal 

rates, and H is a 4 by 1 column matrix for CMG torque. The matrix C has 
no unique inverse. However, the first approximation to the inverse of such a 
matrix is its transpose. The desired control torque, T , based on attitude 

error signals is substituted with a neagative sign for the CMG torque. Approx- 
imating the inverse of the torque matrix by the transpose, the CMG steering 
law is 


< -1 h> C* T 

c 


(103) 


Expanding the transpose steering law, the gimbal rate commands are 
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(104) 


There are no mathematical singularities in the steering law and it 
is easy to implement. The main disadvantages of this type steering law is 
that the desired torque is not produced and the gain through the CMGs depends 
upon the gimbal positions at that instant of time. For example, if the gimbal 
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angles are all zero and the elements of the torque matrix is evaluated at that 
condition, a unit torque command per axis produces a torque of 4s 2 /3 on the 
X-axis and 2c P on both the Y- and Z-axis. More than twice as much torque 
is produced on the X-axis than was commanded, assuming a beta angle of 
53. 1 degrees. However, only about three-fourths of the commanded torque is 
produced on either the Y- or Z-axis. To alleviate this basic deficiency in 
control effectiveness, G. E. [2] has utilized CMG torque feedback in the 
transpose steering law. Since the actual CMG gimbal rate is proportional to 
the torque being produced, CMG tachometers measure the gimbal rates 
which are fed back in a minor loop illustrated in Figure 9. A first-order 

lag filter is installed in the loop to provide added rejection of mechanical 
noise errors. 



figure 9. Transpose steering with torque feedback. 

The modified transpose steering law with CMG torque feedback, k 
becomes 


&1 ' * [ C “ ( T cx - K) * C ” ( T cy - \) * C “ ( T cz - q] /» I 

i 2 = - [c„ (T cx - h x ) + C 22 (T cy - h y ) + C 32 (T cz - h z d /h ( 

r J V . (105) 

‘ - [ c ‘» ( T CX - \ ) + c a ( T ey - \) + C »3 (T cz - hj] /h ( 

' a> ’ - [ C “( T - ' *■*) + ( T cy - l y) + c 3<( t cz - \)} ^ 
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Although the torque feedback provides constant gain for the rate loops, 
the total C MG system still has a variable gain as a function of the CMG gimbal 
angles. However, unless very fine control is required, adequate vehicle 
performance should be maintained. The torque feedback should help provide 
only the control authority required by the attitude error signals. 


G. The BECO H-D i str ibution 

As previously mentioned, equation (22) does not have a unique inverse 
for the CMG gimbal rates in terms of other parameters because the matrix C 
is not square and, therefore, is singular. An exact solution can be obtained 
by finding a constraint equation between the gimbal rates. This equation can 
be used to eliminate one of the gimbal rates from equation (2.2), thus reduc- 
ing the C matrix to a 2 by 2 dimension which has an exact inverse, provided 
it is nonsingular. In lieu of a subsidiary condition between CMGs, BECO [ 2] 
has proposed the following scheme for obtaining a steering law. 

Assume that CMG number 4 is arbitrarily gimbaled at a certain 
constant rate at each instant of time the matrix C is to be inverted. Moreover, 
assume that this rate can be determined just prior to each time step so that 
CMG saturation and ghu conditions are avoided. Under these conditions, 


' v -1 <■ t • (106) 

is used as a fourth equation to augment equations (22) 




Let A be the augmented torque matrix in equation (107) and let A be its 
determinant. Using the elements of C [ equation (31) ] , the determinant of A 
is 

A = (-S/2/h) ^ 2 C 2 /3 Ccy i Co 2 Cct 3 + So , 2 S(oq + aj) 

- Cj3 Co 2 S(oq - « 3 )J . ( io 8 ) 

Assuming A is not zero, the solution of equation (107) is 



where has been substituted for H with a negative sign. 

The solution for CMG gimbal angles numbered 1 , 2 , and 3 [equation 
(109)] depends on both the gimbal rate and position of CMG number 4 . As 
noted by equation (108), the determinant does not depend upon CMG number 
4. The determinant goes to zero under the following conditions: 

1. Q'j = -q ' 3 or «, = 180 - a 3 and oi 2 = 90 degrees. 

2. aq = a ? 3 = 90 degrees. 

3. kj * kg x kg = [ A 4 1 = 0 . 

These conditions can occur without producing ghu. Thus, mathematical 
singularities not coinciding with ghu conditions are introduced by this 
formulation. 

At this point the obvious problem is to suitably determine aq . To 
prevent A from going singular, A must not be zero. It appears that singu- 
larities could be avoided by making the fourth gimbal rate, at each time step, 
inversely proportional to the determinant of A. However, to minimize the 
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total gimbal rates, the fourth gimbal should also be proportional to the root 
sum square of the other gimbal rates. Thus, a , may be determined by 
setting 


Q‘4 


K( 


o i 


O' 2 + O' 3 


V 2 

! ) /A 


( 110 ) 


where K is a constant, A suitable value of K was found to be 0.001. 

There are certain relative orientations of the spin momentum vectors 
of the four CMGs such that no torque can be produced in a particular direc- 
tion, Such orientations will be referred to as gyro hang-up orientations and 
should be avoided if possible to have complete and independent control of the 
spacecraft axes at every instant of time. 

The condition for the gyro hang-up orientation is that at any given 
instant of time, the torque vectors of the four CMGs happen to lie in a plane. 

If at that instant of time a torque normal to that plane is commanded, the 
CMGs cannot contribute any torque in the commanded direction and control 
is lost in that particular direction. Therefore, the only way to avoid gyro 
hang-up orientations is to not allow the torque vectors of the four CMGs to 
go into a plane. The following analysis is conducted to determine the con- 
ditions in terms of the known quantities like a } , a 2 , o . , and a j when 
gyro hang-up orientations may occur and a scheme to avoid these orientations. 

The determinant A vanishes both at ghu orientations and at mathe- 
matical singularities. In principle, if A is positive, the fourth CMG can be 
used to cause A to increase, remain the same, or decrease as slowly as 
possible. Hopefully, this will delay the occurrence of ghu or singularity as 
long as possible. 

To get an explicit expression for a, which will make the derivative 

. 

A positive semidefinite, the following procedure may be adopted: 

1. Differentiate A with respect to time [equation (108)1 to get an 
expression in terms of gimbal angles and gimbal rates. 

2. Substitute the values of the various gimbal rates from equation 
(109) into the above expression for A . The resulting expression is equated 


to zero and solved for o 4 . This expression for a 4 is in terms of gimbal 
angles, spacecraft rates, commanded control torques, and the CMG angular 
momentum resolved into spacecraft body axes. Thus the value of cv 4 is 
completely known at every instant of time. 

It may happen at certain instances that the gimbal rate <a 4 calculated 
from the above procedure will give higher values than the upper limit of one 
degree per second allowed. In such instances, the upper limit value of one 
degree per second for a 4 will be used and in that case the value of A will 
tend to decrease. 

Figure 10 shows an outline of the information flow diagram for this 
CMG control law. The following are important features of the proposed CMG 
control law: 


1. Provides control torques exactly as commanded in all three axes. 

2. No interaxis cross coupling, 

3. Involves only one matrix (4 by 4) inversion. Three of the 
elements of this matrix are identically zero; this further reduces the complex 
matrix inversion computations. 

4. Tends to distribute the momentum between each CMG and, for 
this reason, is referred to as the H-distribution steering law. 

Digital simulations show that the BECO steering law [ using equations 
(109) and ( 1 1 0 ) ] gave performance equaling that of the pseudo inverse; 
however, it is more complicated than the pseudo inverse. For each CMG 
failure, the failure would have to be identified, the BECO steering law 
deleted, and an exact inverse inserted which would depend upon the failed 
CMG. 


SECTION IV. STEERING LAW SUMMARY AND SELECTION 


The problems in selecting a steering law are caused primarily by the 
fact that there are more unknowns than there are equations between the un- 
known variables. For example, the four-skewed CMG configuration baselined 
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Figure 10. H-distribution steering law. 


tor HEAO-C has lour gimbal rate commands as unknowns. However, there 
are only three equations from which to obtain a gimbal rate solution. These 
equations are obtained by equating the three components ol the torque com- 
mand vector to the CMG torque vector. Although there are many solutions, 
an exact solution in the usual sense does not necessarily exist. At this point 
each designer will devise his own scheme for obtaining a solution. What is 
needed, ol course, is a constraint equation between variables so that an 
exact solution can be obtained. For each constraint or assumption that is 
made, a dilterent solution will be obtained for the gimbal rates. Whatever 
the solution, it is referred to as the CMG steering law. Several candidate 
steeling laws were derived and compared on the basis ot their effectiveness 
in producing the desired control torque. Each was derived lor the lour- 
skewed CMG configuration which has been baselined for HEAO-C. and each 
w'as implemented in digital simulations. The candidate steering laws were- 


1. Constant 


gam. 


2 . Maxi mum contribution . 


3. Pseudo inverse. 
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4. The Bendix summation of three-gimbal inverses. 

5. The BECO momentum vector distribution. 

6. The G.E. transpose with torque feedback. 

Given any CMG configuration and steering law, there are certain 
conditions that can produce problems in either the performance of the system 
or in the ability to obtain the commanded torque. As previously defined, the 
CMG torque is related to the gimbal rates by a torque matrix whose columns 
are unit vectors along each CMG torque output axis. Since there are four 
CMGs, there are four columns in the torque matrix. When all four columns 
are coplanar, no torque can be produced perpendicular to that plane. This 
condition has been defined as gyro hang-up. That is, when the commanded 
torque is perpendicular to the total CMG torque vector, it cannot be produced. 
The ultimate in gyro hang-up is CMG saturation. Obviously, if the CMGs 
have produced all the momentum possible in a given direction and the com- 
manded torque vector asks for more in that direction, it is impossible to 
produce the required change in momentum and, subsequently, control is lost. 
For an inertially oriented spacecraft, a component of the gravity gradient 
(s* ) torque is usually offset from its zero reference, producing a momen- 

tum component that tends to saturate the CMGs over a period of time. Before 
the CMGs are saturated, the prevailing condition must be detected and stored 
CMG momentum must be interchanged with that produced by a second source, 
such as RCS or electromagnets. In so doing, the CMG gimbal positions are ’ 
normally returned to their reference position, usually a zero momentum 
state. 


Each steering law was compared on the basis of complexity, accuracy, 
mathematical singularities, failure adaption, performance after failure, and 
growth potential. No attempt was made to weigh the importance of each 
comparative factor and the ratings shown in Table 2 are somewhat subjective. 
However, on the basis of an elaborate digital simulation (Appendices A, B, 
and C) the pseudo inverse consistently gave better performance than the other 
steering laws, especially when various CMG failures were simulated, and 
was relatively easy to implement. The context in which the comparative 
factors were used are as follows: 

1. Complexity — The mathematical manipulation and logic required 
for implementation. 

2. Accuracy — The pointing and jitter performance obtained through 
simulations. 
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TABLE 2. CMG STEERING LAW COMPARISON 



rhe pseudo inverse is recommended lor HEAO-C. 


3. Singularity — Combinations of gimbal angles which can cause 
zeros in the denominator of the steering law, hence, program divergence. 


4. Failure Adaption — The corrective actions that must be taken in 
the event of a CMG failure, in particular, changes in the steering law. 

5. Pertormance After a Failure — Pointing and jitter performance 
with one CMG out. 

6. Growth Potential — Minimum modifications required to accom- 
modate more than four CMGs. 

7. Gyro hang-up — A combination of gimbal angles which prevents 
the desired torque from being produced: (a) cannot transfer attitude error 
through the steering law, and (b) the commanded torque is perpendicular to 
the instantaneous CMG torque vector ( H • T 0) . The ultimate in gyro— 

hang-up is CMG saturation. 

8. Cross Coupling — Not being able to produce a torque about one 
axis without also torquing a transverse axis. 

9. Efficiency — Capability of using the total momentum available 
without gyro hang-up or cross coupling. 

The "constant gain" steering law is derived by assuming small gimbal 
angles and assuming that the CMGs operate similar to scissored pairs. As 
implied by the name, the constant gain steering law contains constants which 
can easily be implemented on an analog computer. It is the simplest of all 
laws but is valid for gimbal angles less than t 90 degrees. It contains no 
singularities. With one CMG failed, the CMG null position must be redefined 
and a new constant gain steering law derived which is valid about the new null 
position. After a failure, the gimbal angles tend to become larger and the 
performance degrades. 

The "maximum contribution" is derived by assuming that each CMG 
operates independently. The gimbal rate of each CMG is commanded to 
produce as much as possible of the desired torque. In complexity, it compares 
favorably with the constant gain steering law. Excellent performance is 
obtained as long as the gimbal angles stay small. However, as the gimbal 
angles become large, gyro hang- up conditions are approached and the per- 
formance is degraded due to cross coupling torques. There are no singular- 
ities in the maximum contribution steering law and no changes are required 
for failure adaption. With only three CMGs operating, the performance is 
degraded due to cross coupling CMG torques. 
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The "pseudo inverse" steering law basically minimizes the norm 
between gimbal rates and is based on the work of R. Penrose [ 8] . A 3 by 3 
matrix inversion is required to get the inverse, along with several matrix 
multiplications. It represents the ultimate accuracy in performance. With 
one CMG failed, the pseudo inverse reduces to an exact inverse without 
program modifications, and its performance is not degraded. Since the 
pseudo inverse tends to produce exactly what is commanded, there are no 
cross coupling torques. However, gyro hang-up conditions can be approached 
as the gimbal angles become large. The pseudo inverse is recommended for 
HEAO-C. 

The "Bendix summation of three gimbal" solutions is much more 
complex than any of the other candidate steering laws. Basically, the CMG 
torque vectors are arranged in combinations of three. There are four possible 
combinations, each having an exact inverse. The solutions to each combina- 
tion are summed to produce the steering law. It is not known what the equiva- 
lent constraint equation would be or what, if any tiling, is being minimized. 

The Bendix law introduces several singularities that are inherent only to their 
scheme. In addition, gyro hang-up conditions can be attained. Singularity 
detection and avoidance techniques are required which compound the com- 
plexity. Failure detection and corrective actions are required. However, 
after a failure only one of the three gimbal inverses would be used, in which 
case the steering law reverts to an exact inverse and the performance 
improves. The Bendix steering law is not recommended for HEAO. 

The "BECO H-distribution" is derived by assuming a constraint 
between the gimbal rates that tend to distribute the CMG momentum to avoid 
gyro hang-up conditions. In addition to a constraint equation, a 4 by 4 
matrix must be inverted, increasing its complexity. It performed well in 
simulations; however, complete reprogramming is required to accommodate 
a CMG failure. With one CMG out, supposedly the exact inverse w'ould be 
used, improving performance. 

The "G.E. transpose with torque feedback" is a variation of the 
maximum contribution steering law'. Basically, each CMG is commanded 
individually with the CMG torque signal being fed back to prevent overcontrol 
and provide stability. It is more complex than the maximum contribution 
with about the same performance capabilities. As the gimbal angles become 
large, gyro hang-up conditions can also be attained. Both the maximum 
contribution and G.E. steering laws offer maximum growth potential. As 
more CMGs are added, the cross coupling between CMGs becomes less and 
the performance improves. The main objection to this type steering law is 
that the system bandwidth varies as a function of the CMG gimbal position. 
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Each steering law was evaluated through digital simulations of vehicle 
pointing performance. Based on factors such as complexity, accuracy, CMG 
failure adaption, performance after failure, mathematical singularities, and 
singularity avoidance, the pseudo inverse CMG steering law is preferred over 
the other steering laws. As a specific comparison, the Bendix steering law 
introduces several singularities that are inherent to their law and make 
necessary a rather complex singularity detection and avoidance scheme. By 
far, the Bendix three-gimbal inverse steering law is more complex than any 
other candidate steering law without a corresponding increase in either 
accuracy or reliability. In contrast, the pseudo inverse is more accurate 
than any of the other laws, adapts to a failure mode without software modifica- 
tion, and performs after a failure (assuming failure identification) without 
any degradation. 

With the magnetic system continuously dumping CMG momentum, the 
candidate CMG steering laws were simulated. In all cases, the stored CMG 
momentum remains near zero; consequently, the gimbal angles stay very 
small. Even under worst environmental torque conditions, the deviation of 
the angles from their null position was only about 3 degrees. As a result, all 
steering laws performed equally well and all produced excellent pointing and 
jitter performance. Since the gimbal angles stay very small, a constant gain 
steering law is adequate for HEAO with magnetic momentum dump. However, 
if rapid slewing is commanded, the gimbal angles became large during the 
maneuver. If HEAO is required to maneuver rapidly (for example, solar 
flare viewing in an antisolar direction), the pseudo inverse or maximum 
contribution steering law is recommended. Moreover, either steering law 
does not require any modification in the event of a CMG failure, providing 
fail operational capability. Assuming one CMG has failed, the magnetic 
system will automatically drive the remaining CMGs to a new null (zero 
momentum) position, without any changes being made in the software. 


As long as the gimbal angles stay small, as with continuous momentum 
dump against the earth's magnetic field by electromagnets, all steering laws 
perform about equally well. As the gimbal angles become large there are 
significant differences in performance. These differences are due to the 
steering laws' ability to cope with singularities, gyro hang-up, and cross 
coupling CMG torques. Based on a comparison of the candidate steering 
laws for HEAO, the pseudo inverse is recommended for I1EAO-C. As an 
alternate, the maximum contribution should be considered, especially if 
more CMGs were added for greater reliability over the two year mission. 
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The HEAO-C control system is described in Appendix A, including 
electromagnet sizing for continuous momentum dump. Euler’s equations for 
HEAO-C with the four-skewed CMG configuration illustrated in Figure 2 are 
derived in Appendix B. The digital simulation of HEAO-C is described in 
Appendix C along with typical performance results using the various CMG 
steering laws. The spacecraft inertia properties, orbital parameters, and 
feedback gains are also given. 
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APPENDIX A. SYSTEM DESCRIPTION 


An alternate actuator system proposed for HEAO-C is composed of the 
four skewed CMGs for maneuvering and attitude hold during normal operations, 
a cold gas RCS for control during orbit adjust stage (OAS) burn and initial 
stabilization, and three orthogonal electromagnets for momentum management 
of the CMGs and for direct control torque in the event of two CMG failures. 

Such a system will be more reliable, weigh less, and provide greater depth of 
failure without degrading system performance than the RCS-CMG system base- 
lined for HEAO-C. For long lifetime missions such as HEAO, a considerable 
amount of RCS fuel is required to dump the accumulated CMG momentum due 
to biased environmental forces. It seems only natural to consider the produc- 
tive use of environmental forces such as gravity gradient or the earth's mag- 
netic field. Previous studies for the Skylab program show that gravity gradi- 
ent can be utilized to dump CMG momentum. However, during dump, the 
spacecraft's pointing requirements must be ignored and the vehicle maneu- 
vered in a specific sequence such that the gravity gradient counteracts the 
accumulated momentum. For HEAO, experiment viewing time takes priority 
and precludes the use of gravity dump except perhaps as a backup measure. 
Currently, the earth’s magnetic field offers the greatest growth potential for 
utilizing the HEAO environment for control purposes, without sacrificing 
experiment viewing time bv imposed maneuvers or restricting the observatory's 
orientation. 


MAGNETIC TORQUER 


Basically, the advantage of utilization of controlled interactions with 
ambient fields is that no fuel need be carried aboard the vehicle for CMG 
momentum dump. However, the use ol electromagnets to react against the 
earth's magnetic field does require additional power to drive the coils. The 
magnetic system proposed for HEAO consists of three electromagnets aligned 
orthogonally with each vehicle control axis which is assumed to be a principal 
axis. When current is passed through the coils, a dipole moment, M, is 
generated^ which reacts with the earth's magnetic field, B, to produce a 
torque, T . The torque produced obeys the vector cross product 
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It is apparent that the torque produced is perpendicular to both the dipole 
moment and the earth’s field. Although limited in magnitude by the electro- 
magnet size, the direction of the dipole can be produced in any direction. At 
any instant of time, the magnitude and direction of the earth's field depends on 
the observatory's orbital position relative to the surface of the earth. The 
vector components ol B would be obtained by onboard magnetometers. The 
magnitude ot M varies as a function of the currents being passed through the 
coils at any time. To maximize the torque produced by a given current, the 
dipole generated should be perpendicular to the earth’s field. Moreover, it 
is apparent that a torque cannot be produced in the direction of B. 

At some instant of time, the desired torque may be aligned with B, in 
which case it cannot be produced. However, these periods are relatively short 
because, as the orbital position of the observatory changes, a corresponding 
change occurs in the direction of the earth’s field. Over any time interval 
during an orbit, the CM Us produce the desired torque required for fine control 
and the magnetic torque, il available, is used to dump the momentum accumu- 
lated in the CMGs. I he magnetic system proposed tor HKAO provides a torque 
proportional to the stored momentum. As such it is a secondary control torque 
and, if it cannot momentarily be produced, the vehicle performance is not 
degraded. 


MAGNETIC CONTROL LAW 


I a- l the CMG momentum he denoted by the vector 

H hi - h j - h k . (A- 2) 

x z v 

I hen, il a proportional system is considered, the magnetic* torque required to 
dump the CMG momenta must be proportional to H but opposite in direction; 
therefore, 

T -K H / a -l 

m m v A -'d 

where K is an arbitrary constant to be determined. Equaling equations (A-l) 
and (A-:s) and taking the vector cross product of B with both sides gives 

B - (-K m iT) B (M B) B-M - (M • B) B . (A-4) 
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The maximum torque for a given magnitude of B and M is obtained when M 
is normal to B, implying that M • B = 0. For this case, equation (A-4) 
can be solved for M to give 


M = -K (B x H)/B 2 
m 


(A-5) 


Equation (A-5) gives the dipole moment required to dump the CMG momentum 
H. In expanded form, the vector components of the required magnetic control 
law for momentum dumping are 


M = 

(-K /B 2 ) (B h 

X 

m y z 

M = 

(-K /B 2 ) (B h 

y 

m z x 

M = 
z 

(-K /B 2 )(B h 
v m x y 


B h ) 

z y 


B h ) 

x z 


B h ) 

y * 



(A-6) 


Of course, the magnetic dipole can be directly related to current and voltage. 
For use on HEAO, the power has been arbitrarily limited to 10 watts per elec- 
tromagnet. Substituting the dipole commands into the torque equation [equa- 
tion (A-l)] produces 

T = (-K /B 2 )(B x H) X B = (-K / B 2 )[ B 2 H - B(B • H)] (A- 7) 
m ' m m 


as the magnetic torque produced to dump the CMG momentum. If H is per- 
pendicular to B, then B • H = 0 and the exact torque needed for momentum 
dump is produced. Consider the other extreme and assume that H is aligned 
with B. In this case, H can be expressed as a constant k times B (H = kB) 
and equation (A-7) becomes zero. That is, no magnetic torque is produced 
when the earth's field is unfavorable for dumping momentum. Only that portion 
of the desired torque which is perpendicular to the earth's field will be pro- 
duced at any given time. However, momentum can be dumped on one axis at 
the expense of increasing momentum on another axis, but the total magnitude 
will always be reduced by the magnetic system. 

In the event of two CMG failures, the magnetic system could be used to 
provide direct torque in addition to dumping CMG momentum. For direct 
torque control, the magnetic torque would be set proportional to the desired 
control torque, T . The desired torque is based upon attitude error signals 
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which have been weighted by appropriate feedback gains. Normally, the CMGs 
would provide this torque through the CMG steering law. The dipole moment 
required for direct torque commands is obtained b.v setting 

K e * c ’ (A-H) 


In a manner similar to that used to obtain the dipole commands for momentum 
dump, the dipole command for direct torque control is 


M K (15 T )/B 2 
e c' 


(A-9) 


Assuming two CMGs have tailed, the dipole commands would be a com- 
bination of that required for momentum dump and direct control. The magnetic 
torque would be set equal to 


T 

m 


-K II + k T 
m e c 


(A- 10) 


The corresponding dipole solution is 


B x (-K II + I< T ) -K 

I'*'* Cl ' 


M = 


m 


K 


B 2 


f (B ■ H) + -g (BxTJ . (A- 11) 


B- 


Hence, the lorm ot the dipole command changes according to the type actuation 
desired. Appropriate values for tne constants are K -0.01 and K =10 

i i 'i e * ' 


Electromagnet Sizing 

Ihe maximum dipole is physically limited by the shape and volume of 
the electromagnet, the number of turns in the coil, current passed through the 
coils, and physical properties ol the materials used. In sizing the electro- 
magnets, low power usage is selected over weight as a design criterion. A 
maximum of 10 watts per coil has been arbitrarily selected as an upper limit 
and the rest of the magnetic system has been sized accordingly to meet the 
required torque and/or momentum dump capability. The magnetic system 
would be installed in the OAS and, for this reason, the length of the electro- 
magnet has been limited to GO inches. 
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Based on simulation results for which the dipoles per axis were limited 
to selected values, it was found that a dipole moment per axis of 0. 2 ft— lb/ 
gauss was adequate to dump the expected secular momentum due to gravity 
gradient torque. However, under worst-case conditions, the magnetic system 
could not dump all the accumulated momentum and the CMGs could saturate in 
about one day. For direct torque control, the magnetic system must produce 
a torque equal to or greater than that of gravity gradient, in which case a 
dipole moment of 0.4 ft-lb/gauss is desirable. With a properly sized magnetic 
system, two out of four CMGs can be failed and still maintain acceptable 
HEAO-C performance. For this reason a value of 0.4 ft-lb/gauss was selected 
as a basis for designing electromagnets for HEAO-C. A candidate electro- 
magnetic torquer design to meet the above specifications is shown in Table A-l. 
In this case, an AEM 4750 Core was assumed to be utilized. 


TABLE A-l. ELECTROMAGNET TORQUER DESIGN DATA 



Total 

( 1 torquer) 

T otal 

( 3 torquers) 

Weight, lbm 

110 

330 

Max Power, W 

10 

30 

Outside Diameter, in. 

2.45 


Core Diameter, in. 

2. 11 


Core Volume, in. 3 

209.5 

628.5 

Max Magnetic Moment, 
amp-turn-m 2 

5440 

9422 

Torque Produced in a 0.35 
gauss field, ft-lb 

0. 14 

0.24 

Flux Density, gauss 

12 000 


Field Intensity, oersted 

20 


Core Material 

AEM 4750 


Winding Material 

Aluminum 
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As simulation data were obtained using coils for CMG momentum 
dump, an improvement in observatory performance was noted. A linear 
analysis of the HEAO-C equations of motion proved that magnetic momentum 
dump inti oduceci the integral of attitude error through the control loop, improv- 
ing both pointing and jitter performance of HEAO-C. Moreover, since momen- 
tum is continuously dumped, the CMG gimbal angles stay very small, permit- 
ting the use of a constant gain steering law. Since the coil commands are 
based on the CMG momentum state, when one CMG is failed, the remaining 
CMGs are automatically driven to a new' null (zero momentum) state without 
reprogramming or software modification. The following are some of the 
advantages of using electromagnets lor continuous momentum management: 

1. No fuel or RCS required for momentum dump. 

2. Lifetime not limited by expendables. 

■ Saturation detection not required. 

4. Very small gimbal angles permit use of a constant gain steering- 

law. 

5. One CMG fail operational capability. 

fi. Operation with two CMGs failed is possible. 

7. Small size CMGs (50 ft-lb-see each) could be used. 

S. Improved pointing performance". 

there are, however, some possible disadvantages in the use of electro- 
magnets. These are possible magnetic contamination, which may require that 
certain components such as photomultiplier tubes be shielded, or power usage, 
which would be limited to ;]() watts for the three coils. Overall, the system” 
should be more reliable than one using KCS dump. 

SYSTEM COMPONENTS 


1 he altitude sensing and control system components are shown in Fig- 
ure A-l. I he sensors provide attitude information which is processed in a 
central computer to generate attitude error signals. The control laws and 
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Figure A-l. HEAO-C attitude senping and control system. 

algorithms are generated within the computer to drive the actuator systems. 
The actuators (CMGs, RCS, and electromagnets) provide the torque required 
for attitude hold against environmental disturbances and for maneuvering. 

Most of the components are redundant so that fail operational capability exists 
for most subsystem failures. As illustrated, there are three coarse sun sen- 
sors, two digital sun sensors, six rate gyros, four fixed head star sensors, 
and two 3-axis magnetometers. There are three input-output processors and 
computers with only one operational at any time. The RCS is completely dual 
redundant with only half the system normally in use. There are three orthog- 
onal bar torquers in the magnetic system, each with separate drive electronics 
and all are normally operating. All four skewed CMGs are also normally 
operating. However, any bar torquer or CMG can be failed without degrading 
system performance. If two CMGs fail, control can be maintained if the mag- 
netic system is used to provide direct torque as well as CMG momentum 
management. 
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The interaction between the attitude sensing and control system and the 
HE AO spacecraft is shown in Figure A -2. Ground commands for pointing to 
various target source's are used to program and/or drive the guidance, navi- 
gation, and sequencing logic in the control computer. Based on the attitude 
commands and the present vehicle attitude from the onboard sensors, the 
vehicle control law forms an attitude error signal, usually denoted as the 
commanded torque, T . Then, the commanded torque must be processed 

by an appropriate actuator steering law to obtain signals to drive the actuators. 

The actuators (CMGs, RCS, or magnetic torquers) produce a control torque 

on the vehicle that counteracts the disturbance torques, T , and forces the 

d 

attitude* errors to zero. Based on the ("MG accumulated momentum status 
and the earth’s magnetic field, the magnetic coils are energized to produce 
a dipole moment which reacts with the earth’s magnetic field to produce a 
magnetic torque, T . The* magnetic torque is normally used continuously to 

keep the CMG momentum near zero, thus, preventing C MG saturation. With 
the magnetic* system, tin* RCS is not used in normal operational modes after 
the orbit has been established. Either with or without the magnetic system, 
however, an RCS is still required for initial stabilization, control during 
orbital adjust stage burn, initial solar acquisition, establishing the first 
celestial reierenee, and attitude* hold during checkout before and while the 
CMGs are activated. While establishing the orbit and during initial checkout 
phases, the RCS is used continuously as required. After the CMGs are oper- 
ational, the* RCS is only used periodically for CMG momentum dump or 
emergencies. Normal maneuvering and attitude control is done by gimbaling 
the CMGs. 
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Figure A-2. Simplified attitude control system block diagram. 




















APPENDIX B. EULER'S EQUATIONS FOR HEAO-C 


’l hr dvtutmie equations which govern l he rotational motion of a rigid 
II KAO with four skewed single gimbal CMGs are obtained by equating the time 
derivative ol the total system angular momentum to the sum of the applied 
torques: 


H V T(applied) - T - T 
t m g 


(B-l) 


The applied torques are due to environmental forces such as gravilv, aerodv- 
namie and electromagnet interactions, and to onboard actuators such as reac- 
tion jets which expel mass from the spacecraft. Simulations have shown that 
for higher orbits, the dominant environmental torque is gravity gradient, T 

( T 

and the others can he ignored during preliminary design studies. Normally, 
torque due to electromagnet interaction with the earth’s ambient field is very 
small. But, since* an electromagnet system is proposed for continuous momen- 
tum management of the CMGs, the magnetic torque given by the equation 


must be included as an applied torque. Although a reaction jet control system 
will he < e m 1 is tun eonudered as a pari of the huler equations 

since the nasn ubjoclnv i.-> to show the performance of the- GIG system along 
with bn r loi < |> ; < , 


The total s\ stem angular momentum is composed of two parts: that due 
to the npaeecralt motion and that due to the CMGs. The spacecraft angular 
momentum li^ i s tne product oi u.s inertia tensor, 1, and angular rale, 

\ 1 ’ • (B-.J) 


For slurb purposes, the products of inertia are assumed to lie zero, in which 
east* 


n 


\ r 


z z 


(B-4) 
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is the vector form of the spacecraft momentum relative to its reference coor- 
dinates. The CMG momentum has previously been derived in the spacecraft 
reference coordinates as: 

H = h i + h j + h k , (B-5) 

c x r y r z r 

whose components are given by equation (28) . The total vehicle momentum is 

H = H + H . (B-6) 

t v c 


The time derivative of any vector relative to inertial space is equal to 
its derivative relative to its reference coordinates plus the cross product of 
the angular motion of the reference coordinates relative to inertial space and 
the vector. The angular motion of the spacecraft body reference coordinates 
is given by 

to = cui+cuj+cuk . . (B-7) 

x r y r z r 


Selecting body reference coordinates as a basis in which to perform the vector 
operations, equation (B-l) becomes 



inertial 



reference 


+ 


cu x H 


T + T 
m g 


or 


* * * 

H + oSxH + H + cu X H = T + T 
v v c c m g 


(B-8) 


(B-9) 


Performing the indicated vector operations produces the Euler equa- 
tions for HEAO-C. The vector components of equation (B-9) produce three 
equations which govern the rotational motion: 
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As derived in the main text, equations (28) and (30), 


and 
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(B-12) 


I 

1' •'(! - i ) ii it \ 
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(B- 13) 


ii , U , ; , hd u ni r ; a. , ( -i v,\ )u;u. uts ul 1 he local vertical unit vector pro- 

: a, : . a ► , (, > « :,'••• ri ioi'> If U‘s and f is the orbital angular rale. 

o 

a i b I < * n»i .in- a; eq; ;utii>n (B-2), are 

1 Aii > 

m.\ \ 

1 .VI n - A I i- ; , (ii-il) 

i > ■ S / ' v / 

1 A! i, \i [; 1 

Ilf . V V \ / 

\1 , AS . M are die e!( e( i oiivuii'iel dipole components and li , I> , 

\ \ /. \ v 

arc the earth’s magnetic Field components, all expressed in Inxlv 



coordinates. 



The spacecraft is related to solar inertial coordinates by a Euler angle 

sequence. The solar coordinates are defined by having the X -axis pointing 

s 

to the sun, the Z -axis perpendicular to the ecliptic plane directed northward 
and the Y g -axis completing a right-hand triad in the ecliptic plane. The trans- 
formation between coordinates is derived by first rotating about the X -axis 

S 

by the angle 0 , then rotating about the once transformed Y’-axis bv the 

s 

angle 9 and, finally, about the twice transformed Z”-axis by the angle ip . In 

s 

vector matrix form, 

^r ^rs^s ‘ (B-15) 

The elements of B (a 1, 2, 3 rotation order) are 


Bn 

= CO C 0 


B 12 

= C0 S ip + 

S 0 SO Cip 

B13 

= S0 Sip - 

C0 SO op 

B 21 

= -CO sip 


B 22 

= C0 C ip - 

So SO Sip 

B 23 

••= S0 C Ip - 

C0 SO Sip 

B 31 

= SO 


B 32 

= -S0 C6 


B 33 

= C0 CO 



/ 


(B-16) 


The body angular xates, equation (B— 7), are related to the Kuier angular rates 
by the kinematic relations 

co = (pCOCip + 9 Sip j 

x i 

Uy = 9Cip - 0C0S ip > . (B- 17) 

= 0 + 0 SO J 
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The inverse relations are 


U' _ r - - x S;:)/('n | 

"\ tV * • L 'x S - / • (B- 18) 


Addition.)! t nonimal: s and t ransformations between coordinates are required 
to generate iiit earth magnetic lield and the local vertical vector components. 

I heso arc conlained in llu 1 1 KAO- A summary report NASA TMX-53976 | 1| and 
ilie simulation tic limiton documeiu prepared by Computer Sciences Corporation, 
Huntsville, Alabama |9|. Also, additional iniormation on the transformations 
a ad gravity lorque :■ »*«.• contained in NASA TMX-58829 | 10). 

1 or stwdv purposes, the HKAO body reference coordinates are defined 

by X v be b i - | erpcndieida r to the hard-mounted solar panels, Y which is 

r 

along the axis ol minimum inertia (long axis), and Z which completes a 

tight-hand triad. When the 1 uler angles are zero, the body reference and 
solar coordinates are aligned. Vehicle pointing performance is measured by 

the solar oliset angle, 6 , the target pointing error, 6 , and the spacecraft 

s P 

angola i * rat . . -■ . 


; . ingle i s given in 

i , 


\. ‘'iW (B n ) * (-• • „<-) 


ltd the i:i; 


>. ; ;o: m i ng error I )\ 




(H- 19) 


(B-20) 


W lie IX ,) * 

point inn iim 
spaeeo ra i ! 


■ - ;!iM ^ ! h< rom man* It. d I uler angles re([uired ior target 

i,,n - |> icccrali axis loan experiment target source. The total 

I ibilit' , commonly referred to as jitter, is given bv 


(H-21) 


bir roll 


,lU ’ iinc-oi-sigat, long observatorv axis, is given bv 


6 I o o 

r c 


(B-22) 


77 



During normal operations, the HEAO-C performance requirements are 


5 < 37 degrees 

s 

6^1 arc minute 
P 

5 < 5 arc minutes 

r 

u) < 1 arc second/ second 



(B-23) 
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APPENDIX C. PERFORMANCE SIMULATIONS 


INTRODUCTION 


A fliiiii.il computer program was written to simulate the dynamic behav- 
ior ol l lie IIKAO-C spacecraft in a circular earth orbit. The program is basic- 
all v a modification of the program [9| used for IIKAO-A with the addition of 
('MCI dynamics and steering law. The equations which were programmed were 
I.uier’s equations lor rotational motion about the principal observatory body 
a.\es, buler s kinematieal relations which relate the observatory body princi- 
pal axes to the solar reference, transformational matrices which relate the 
environmental forces to the observatory body axes, control logic which relates 
t ie llbAO spacecraft s attitude errors and rates through appropriate feedback 
gams to applied torques about the observatory body axes, a spherical harmonic 
expansion ol the earth’s magnetic field, CMC dynamics, several selected 
steering laws, and magnetic control torque logic for CMG momentum manage- 

mtMil. 


' <n ‘ ul(l ' ,( l icalism, the program includeo all natural movements wnich 
could aik it the spueeeialt s attitude motion. These movements include the 
earth’s revolution about the suit (1 deg/da\ ), regression of the ascending line 
oi orbital node ((> deg/da.v ) and the earth's rotation (:!«() deg/day) . The environ- 
mental forces acting on the spacccrait were gravity gradient, aerodynamic, 


Urn dominant cm i ronim-nla I torque ‘is that caused by g.g. effects. To 
simulate the gravity torque, the vector components of the local radius vector 
must be projected into body coordinates. The effects of orbital position, incli- 
nation, orbital regression, time of year, and position of the ascending line of 
nodes were considered in deriving the required transformational matrices. 

1 he magnitude ol the g.g. torque also depends upon the vehicle inertia prop- 
erties, m particular the difference between the inertia values. Both spacecraft 
inertia values and orbital conditions were selected such that the g.g. torque 
attained its maximum value. 


. 1 1UUl l> ( f illustrates a sjinpJitied block diagram of the HKAO-C digital 

simulation. I he simulation contains lour nonlinearities: ( 1 ) limits on the 
(’.MG ginil.nl rales; (2) limits on the CMC. gimbal positions, which are not 
" ,!OU ''’ l [ ' ) 1 ;u-1 ' l ,<>: ibon I cedi. a ek channel (not shown); and (-J) limits 

ou Hie c leet roniagnet dipole- moments. During the simulation runs, the limits 
■ n tue ( MG gimbal positions were ueiet -<l sim e tne CMGs were assumed to 
have unlimited angular positions. 
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Figure C-l. Simulation block diagram 
















Parameters and Cases Simulated 


i h, j : ! ■ *: 1 1 r j ■; u ;--. utilized in Hit' simulation lor base runs are listed 

i z i ! ) s sing-iPy lon»' (pointing) axis inertia 

I S2 1 15) slug -fix intermediate inertia axis 

! X2 r.. ; n .duy-fl % maximum (sunward) inertia axis 

•\A 2i’n ii.iiii., orbital altitude 

i 2s. f> dew, <; rbi tal incd i nat ion 

\ 2 V ' > :!(•!:, winter sok-tiee lime oi year 

ii iso deg. mi)! , ni!i:> terminator line ot nodes 

il j > t|i. w, position of ( . reemvi eh meridian relative to Aries 

ii 2du 1'i-lb- see, momentum per (\M(J 

! o : . t .■.!'•■ ! dipole limit 
■ j. ad \ - a . : i s no -alien P vdl >a ek ga I n 
I. sup/. Y axis position teedbark gain 

i - ■■ ! t n > / peso Pm P edbaek ga in 

i\ , d I 12, X-axi.^ rale leedbaek gain 

j l no, \ -axis ran leedbaek l; a i ti 
- a> /.-axis rate leedbaek gain 

i. id It - ib rad, X-axis position leedbaek limit 

i. ■ i 7 s It b. i.id, Y-axis position leedbaek limit 
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L =481 ft- lb/ rad, Z-axis position feedback limit 
zm 

=1 deg/sec, CMG gimbal rate limit 

K =0. 01, magnet dipole gain for momentum dumping 
m 

/3 = 53. 1 deg, CMG skew angle 

K =1.0, magnet dipole gain for direct torque if two CMGs have 
failed 

More than 100 cases were simulated for various HEAO-C configurations, orbi- 
tal conditions, and control gains. Utilizing the base run parameters, with the 
g.g. torque near its maximum value, six CMG steering laws were evaluated. 

On the basis of this evaluation, the pseudo inverse steering law is recommended 
for HEAO-C. Most of the simulation results shown are based on the pseudo 
inverse CMG steering law and all on one set of the base run parameters. Dur- 
ing the study period, several configurations were evaluated. These were a 
HEAO configuration with distributed subsystems, a configuration with a sub- 
systems module, and several configurations which had the orbital adjust stage 
(OAS) attached. The inertia values shown above represent a growth version 
of the HEAO with distributed subsystems and with the OAS attached. 

The feedback control gains were selected to give a damping ratio of 
0.7 and a natural frequency of 0.314. The corresponding time period is 20 
sec and the system time constant is 10 sec/ rad. As shown in Reference 11, 
the pointing error is proportional to the disturbance torque magnitude and 
inversely proportional to the position feedback gain. In essence, the greater 
the required pointing accuracy, the higher the feedback gains. Introduction of 
integral position feedback permits the use of lower position and rate gains. 
However, HEAO-C pointing performance can be obtained with only position and 
rate feedback terms in the attitude error signal. As an alternate, electro- 
magnets are used to dump the accumulated momentum against the earth’s mag- 
netic field. In this case, the magnetic loop also provides integral feedback of 
the attitude error signal which improves vehicle pointing performance as an 
added bonus. The objective of the attitude error signal is to formulate some 
desired vehicle torque command. 

The objective of the CMG steering law is to convert the torque com- 
mands into CMG gimbal rate commands. When the CMG gimbal positions are 
moved, a corresponding change occurs in the CMG momentum vector, hence, 
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producing a control torque which acts on the vehicle to counteract the disturb- 
ing torque due to the environment. Idea II v, the steeling law should make the 

( M 'is produce exactly the torque which is commanded. 

AH the simulation data ha\o been plotted as a Junction of orbit time as 
measured hum the ascending line of nodes, which tor all cases shown is at 
the morning terminator. At a circular altitude of 270 nautical miles, the 
orbital ' Mil' i- 1. 10702 10“' rad/see with a period of 50S(j see. The orbital 

lime is sealed in 1000-see increments on all graphs. The time step used in 
the digital simulation was 0.5 see; that is, at each half-second of orbit time, 
the equations in the* simulation were solved. The printout was on each 50-see 
internal. Although not important for altitude control purposes, the minimum 
period oi orbital darkness is alxuit 20 minutes pet' orbit and the maximum 
a bon 1 50 minutes. There art 1 about 15.2 orbits pea' day. 

i'li t vehicle peiTormanca is measured In three Kuier angles: <p about 

the simlinc {X -axis), 0 about the once transformed solar Y -axis, and i! 
s s 

■‘•bout the Jwict t ranslonnod solar X -axis, f or small angular deviations the 

s 

!' •sder angles ■;> , ■ , and f correspond to rotations about the spacecraft Z , 

X , and v i axes respectively, The pointing performance is measured by the 
angit b« * l ween the long spa eee ral l axis (X ) and the desired target and is given 

bv 



via t ■ ih = ub;-k j’ipt c indicate.*! dr < -oni : na n< led relations required to point to 
the sehoJed larg'd. file roll error is about tilt' target 1 ine-ol -sight (experi- 
ment a *.:i s ) : 


The poi ni i nr stability '{jitter) is gi\'en i:y 

( U' " X ' } 

.V 

wlievr i he bod\ angular rate. To mo 1 ! II 1 . A ( ) } >e rfo r manee requirements, 

i he e i - 1 or.. must lx hms than 1 arc min poinling, 5 arc min roll, and 1 arc sec/ 

' ec ) d : % 1 ; ■ , 
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Response Capabilities 


The maximum vehicle maneuver capabilities are determined by both 
the CMG system torque and available momentum for making the maneuver. 
Assuming 250 ft- lb-sec per CMG with the gimbal rate limited to 1 deg/sec, 
each CMG has a torque capability of 4.36 ft-lb. Utilizing a four-skewed CMG 
configuration with a skew angle of 53. 1 degrees, the maximum torque per 
vehicle axis is about 14 ft-lb on the sun pointing axis and 5 ft-lb on the two 
transverse axes, as shown in Table C-l. The maximum rotational rate 
' imparted to the vehicle is 2. 66 rpm about the axis of minimum inertia 
(X-axis). On the sun pointing (Z y ) axis, the momentum capacity of 800 ft— lb- 

sec allows the CMGs to spin up the vehicle to 0. 109 rpm. Several runs were 
made with HEAO-C operating in a spinup mode similar to that required for 
HEAO-A. In all cases, HEAO-A performance requirements [ 1] were attained 
with a wide margin. Maximum CMG torque on the Z^-axis permits rapid move- 
ment about the sunline for initial spinup or for small maneuvers about that 
axis as in normal HEAO-C operation. In case of unusual maneuvers, as in 

solar flare viewing, the spacecraft could be rotated about the Y -axis bv 90 

v J 

degrees. If a time optimal maneuvering command were issued, the rotation 
for this antisolar viewing would take about three minutes plus another three 
minutes estimated for settling out time to attain HEAO-C pointing require- 
ments. This gives a total time ot about six minutes. However, a time opti- 
mal maneuver scheme was not incorporated into the HE AO simulation and this 
is an area for future study. Normal maneuvers of 90 degrees or less were 


TABLE C-l. VEHICLE RESPONSE CAPABILITIES 


Axes 

CMG Momentum 
(ft- lb- sec) a 

CMG Torque 
(ft-lb) b 

Max Turning 
Rate (rpm) c 

Time for 180 deg 
Turn (min)^ 

X 

V 

800 

5. 24 

2.66 

2. 18 ! 

Y 

V 

800 

5. 24 

0. 110 

6.54 

Z v 

800 

13. 95 

0. 109 

6.59 


a. Four 250 ft-lb-sec CMGs skewed at 53. 1 degrees. 


b. Max torque with gimbals at null; 1 deg/sec gimbal rates. 

c. At CMG saturation. 

d. Time optimal maneuver, settling out time not included. 
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satisfrmtori I v conducted through t bc^ IJKAO simulation by position commands 
oni\ * baeh position n edutek channel was limited to 20 deg'/ min equivalent 
spurt eratt rotation \n ‘ • ; i s . Tin- tinu* repin i a*d to go from one X-ray source 
to another ! short tot mail angle mancuviT:. and the CMGs do not have time 
to i each a usdioj. Normal pointing mancmers can be made by position 
commands only, without utilizing a more complicated time optimal scheme. 

Orbit Environmental Conditions 

All mm an mem ea^es were run with the orbital parameters as variables, 
the si mu In l na r« : nits m be shown are all lor onlv one typo orbit so that each 
sterna ! r- r;s: : -oiojrrvrl under the same en\ i ronmental conditions. The 

orbit uivl v ' ,M • . . ■ a odious re lative to the solar reference coordinates 
are shoe:, e ; ‘ w . 'll, ■ run is al its winter solstice position (\ 270 

degrees), the m - : hi.: line r no< km is at tla* morning terminator (S2 180 

degrees), and ;lr ms a >> (nut ssa ' ai) i j ; t - its non h magnetic pole tilled toward 
the- r • i i n tV : ... >1- < ■. i . » a. spam i i ad i> oriiaUed with its Z -axis 

' u 1 V 

pointing to u . a;. r : » u m 1 i ■ e : rm ! ad northward perpendicular to the 

ecliptic pi a n> • , and it-: i ad; e- anpO 1 i ns a triad in the ecliptic plane. With 

this Pane urbd ad \da. h- h dm(m son, dm gi avil\ gradient disturbance torque 
ia a, 'nuauum component appears on 

th* .. < ■ : • • . , ■ • i a ; a in a • .. ■ si mi, I- is the gieatcsl potential tor 

gj .. . ’ i. v .. a < - sm . 1 : .....mddi<m will sat uralv the CMGs 

qu i C 1 i i h . i . l s i ■ \ ; i i i < f i U . \ J e p 1 i. S w i i i i a ■ \ \ , 1' : > l 1‘ilSu lor aliitude hold. 

d a- a v -a d lisi li < ■ ■ I he a. a. u.pi M vector (components and maglli- 


tude rein:. i • i.o i a = -o - • , 

jam. i . a . ■ i a i i 

v , s in fl-ll). 'The -c-omponent is 

biaau ' W a 1 . ' : : ! O' - : a - c . 

however, ea , cause e i , ! 

s • : S ;i aiii- • l 

hast '.(1 1 O 1 1. {Ul ; ( ■ 

• il Id. ( H ht r orbital conditions, 

t a r mi either tht‘ X - or Z -axis, 

or to be pro-- p : u. = \ 

a. ;H CM; a-’ Si mi, 

1; •oisclx. Both ih<- X and Z axes 

V V 

have t i tv ■■ - ! ■ i m m- a • i 

>m , ( inc to Veil 1 

>d 1 ,nm at v, the X -component is 

near zero. The nmgmba 

pointing e r i < >r is oi r. cl 1 ■ 

!c Ul 1,1 j - ' s.g, ! 
propo id tuiiai l 

orque is a.duLit 0.1(5 li-lb. 1 lie \ehicle 1 
o i i . . With a control system that 

performs vail, me puinti 
The i »bse r\ al o i v muc u ni 
Figure; C-;; .ire .dr > w n in 

ing errin' takes 
urn com] h au.aU>- 
I i g » 1 1 e i - i , s 

an die shape of t ho disturbance torque, 
, i i i responding to the torque shown in 
M. ia -axis s\minetry both the g. g. 

torque and momentum a r 

i • nea r /- re on 

that axis. Idle cvelic Z -axis torque 
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produces biased but cyclic momentum. However, the biased Y^-axis torque 

produces a linear buildup of momentum with time that will eventually saturate 
the CMGs. This momentum buildup over a period of time is denoted as secular 
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Figure C-2. Solar reference, standard orbital position 
and spacecraft orientation. 








momentum and must be dumped to prevent CMG saturation. As shown the 

momentuJ11 is 405 ft-lb-sec per orbit and, as previously shown by the 
CMG momentum envelope, lour 250 ft-lb-sec CMGs produce a maximum cap- 
ability of about 800 ft-lb-sec in the Y^-axis direction. Therefore, if momen- 
tum is not dumped, the CMGs will saturate in about two orbits. Computer runs 
verified the CMG saturation time. If, however, the CMGs' momentum size 
were doubled to a value of 500 ft-lb-sec per CMG, then, under worst-case 
conditions, about four orbits are required to reach saturation. One general 
problem area for future study research is to investigate techniques for dumping 
only the secular momentum without also dumping the cyclic portions. Starting 
with zero initial conditions, the CMG momentum components are identical to 
the g. g. momentum up to saturation, at which time the two diverge. 


The Earth's Magnetic Field 

The earth's magnetic tield is generated by a subroutine, "B-Field. " 
For proper operation, the spacecraft's latitude and longitude relative to the 
Greenwich meridian must be input at each time of calculation, along with the 
altitude and the time of year. The B-Field outputs are the vector components 
and magnitude ol the earth's magnetic field relative to a geocentric coordinate 
system (eastward, southward, and outward directions) . The parameters A, 
S2 , and , and the orbital position are used to calculate the latitude and 

longitude. The B-Field outputs are operated upon by appropriate transforma- 
tions to obtain the field components in solar and in spacecraft body reference 
coordinates. Essentially, the B-Field 2 and its subroutines perform the 

functions of a magnetometer for measuring the earth's magnetic field in bodv 
coordinates. 


Figure C-5 gives the components of the earth's magnetic field for the 
standard type orbit predominantly used during the simulation. The year 1974 
was used as a reference for calculations. The southward component is always 
negative since the field dipole is directed from south to north. Near the 
magnetic equator at orbital times 0, 2800, and 5700 seconds both the eastward 
and outward components are near zero. Since the orbit starts at the morning 
terminator and the north magnetic pole is tilted toward the sun, the maximum 
outward component, B f , is attained at an orbital time corresponding to 


2. The B-Field, with its spherical harmonic coefficients, digital program 
can be obtained from the National Space Science Data Center, Goddard Space 
Flight Center , Code 601, Greenbelt, Maryland 20771. 
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B-FIELD (GAUSS) 



Figure C-5. B-Field in geocentric coordinates versus orbital time < 10 3 sec) . 

about one-fourth ( 1400 sec) and three-fourths (4400 sec) orbit. Since the 
spacecraft's latitude (Fig. C-0) only attains 28.5 degrees (equal to its orbital 
inclination) the southward component never goes to zero. However, the east- 
ward component is always near zero. The spacecraft longitude did not attain 
a value of 180 degrees during the orbit because of the earth's rotation of 
about 23 degrees per orbit. Both the spacecraft's longitude and latitude are 
shown in Figure C-0 with the corresponding B-Field outputs shown in Figure 
C-5 in gauss units. 


The B-Field components shown in Figure C-5 are projected into solar 

coordinates and from there to observatory body reference coordinates by the 

Fuler angle sequence shown in Appendix B. The earth's magnetic Held com - 

ponents, as would be measured by an onboard magnetometer, are shown in 

Figure C-7. The sunward component B is, in this case, alwavs negative 

z -a 

and attains a magnitude of -0.39 gauss, the northward (perpendicular Lo the 
ecliptic) component, B , is slightly biased positively with a maximum value 


of 0.26 gauss, and the B^ component in the ecliptic plane perpendicular to 

the sunline is near-cyclic. The maximum field magnitude is 0.43 gauss at 
4400 sec and the minimum is 0.24 gauss at 2700 sec. The minimum occurs 
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EARTH'S MAGNETIC FIELD (GAUSS) LAT. AND LONG. (DEG ) 


Rli] 



Figure C-6. Spacecraft longitude and latitude versus 
orbital time ( 10 3 sec) . 



Figure C-7. Earth's magnetic field components in vehicle 
coordinates versus orbital time (10 3 sec). 


m 


when the spacecraft crosses the equatorial plane which is also near the 
magnetic and equatorial plane line of nodes. The maxima occur when the 
latitude is greatest, at about the one-lourth and three- fourths points of the 
orbit. The components given are for a standard reference. However, by 
reorienting the vehicle, completely different values could be obtained. But 
the total magnitude curve, B , would remain the same regardless of 
orientation because it depends only upon altitude and position of the center of 
mass in orbit. 


CMG STEERING LAW SIMULATION 
The Constant Gain Steering Law 

1 he cor slant gain steering law was the one first implemented because 
oi its simplicity. The CMG gimbal rate commands are 
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for a skew angle of 53. 1 degrees and 250 ft-lb-sec CMGs. Because the gains 

K and K have been derived to linearize the system when operating about 
A B 

the CMG null position, pointing performance is expected to deteriorate when 
the CMG gimbal angles get large. Figure C-8 illustrates typical performance 
obtained by using the constant gain steering law. As long as the gimbal angles 
are small, the pointing error is about equal to that obtained with the pseudo 
inverse. On the first half-orbit the peak error is less than 0. 07 arc min and 
the gimbal angles (Fig. C-9) are less than 45 degrees. As the gimbal angles 
become larger during the second half of the orbit, a corresponding increase 
occurs in pointing performance. At about 4400 seconds, gimbal angle number 
3 attains a value of 90 degrees at which time control is lost and the system 
diverges. At this time, the CMGs have only accumulated 338 ft-lb-sec 
momentum (Fig. C-10) . Using a constant gain steering law, less than one- 
half of the total CMG momentum can be used for control purposes. For the 
HEAO inertias and CMG size, momentum would have to be dumped each 
half-orbit to prevent the gimbal angles from becoming too large. By increas- 
ing the momentum per CMG to 500 ft-lb-sec, excellent control over a full 
orbit under worst-case conditions could be maintained without dumping CMG 
momentum. 



Figure C-8. Constant gain pointing performance. 
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Using three orthogonal electromagnets to generate a magnet dipole 
which reacts against the earth's magnetic field, continuous CMG momentum 
dumping is provided. As shown in Figure C-ll, the pointing performance 
improves considerably by using continuous momentum dumping by magnetics. 
The peak error is 0. 022 arc min at 4500 seconds. The gimbal angles (not 
shown) deviate less than 4 degrees from their null position. As long as the 
gimbal angles are less than 15 degrees, the constant gain steering law 
performs as well as the most complicated laws or the pseudo inverse. The 
accumulated CMG momentum is shown in Figure C-10 for the constant gain 
steering law, both with and without continuous momentum dump. Without 
magnetics, the momentum value is 338 ft-lb-sec at 4400 seconds after which 
time control is lost, the vehicle rotates violently, and the momentum is 
absorbed by increased vehicle angular rates. With magnetic dumping, the 
accumulated CMG momentum is less than 30 ft-lb-sec, indicating that very 
small CMGs could have been used for pointing control. 
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Figure C-ll. Constant gain pointing error with magnetic momentum dump. 


Maximum Contribution 

The maximum contribution (MC) steering law has been derived such 
that each gimbal rate is commanded independently of the other gimbal rates 
or positions. The steering law is 
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and 


K -0.UU554S15 

are constants that have been tlerived to provide unity gain through the CMG 
loop. 


1 ' tne MC steering law. each CMG gimbal rate is commanded 
Gis d *■ *:lon?* 'V‘ )'(■ •'•onl.ro! ] ing) to contribute as much as possible to the 
commanded torque. Simulations indicate that when the gimbal angles get 
large, certain combinations of angles can momentarily prevent the desired 
torqm h-mi being produced without also producing unwanted torque compo- 
nents that will disturb the pointing accuracy. This is indicated in Figure C-12 
by run number MC-68: All the CMGs are operating and no momentum is 
being « in liiped. At 16UU seconds, the pointing error is 1.26 arc min winch 
exceeds IIKAO-C pointing specifications. The roll error attains a peak 
amplitude oi 4.45 arc min. However, by using higher feedback gains, the 
error can be reduced to an acceptable value (shown by a later rim not included 
h rc h ) . Ait* r (he peak m’or. acceptable performance was maintained 
>' 1}f • ; b In r :H irat !oi« near the t •. o * u I >i t lime poi ni. Near 1 1 • e peak a l ti tilde 
c* r< m . tin- gimbal angle:-. is > w n in figure C 1:;. are changing rapidly. 

\!lcr ni.e mini. the third gimbal angle attains a value of 142 degrees. The 
j u..je- si <* deviate* cuimulcrably Ji'oin their mill position. The CAIG 
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RUN MC-68 



ORBIT TIME (10 3 SEC> 


gimbal rates are all hard limited at 1 deg/ sec (0.017453 rad/sec) . However, 
the limit is not approached during normal attitude hold modes. The gimbal 
rates shown in Figure C-14 correspond to the angles previously shown in 
Figure C-13. These rates are typical for all steering laws (except the trans- 
pose with torque feedback) . hi this particular run, a peak value of 0.055 deg/ 
min, occurs at 4500 seconds. 



Figure C-14. CA1G gimbal rates with maximum contribution steering. 

CMC number 3 was failed by setting its momentum to zero. However, 
no change was made in computing the gimbal rate command for the remaining 
CAIGs. They were initialized to a new zero momentum state by setting 
o? j = : 0, Q?o = 56.4, and a j = -56.4 degrees. The resulting pointing per- 
formance is shown in Figure C-15. The three CMG system hit a gyro hang-up 
condition at 4100 seconds after which pointing control was lost. Before satura- 
tion, both the gravity gradient and stored CMG momentum are identical but 
differ once the CMGs have saturated. As shown in Figure C-16, about 300 
ft-lb-sec have been stored in the CMG system when control was lost. The 
momentum envelope for the three CMGs shown in Figure 7 has an indenture 
centered around the gimbal axes of each CMG, at which point the maximum 
momentum is only about 368 ft-lb-sec along a CMG gimbal axis. The CMGs 
have been driven into ghu near this minimum momentum state. 
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RUN MC-70 
CMG #3 OUT 



ORBIT TIME (10 3 SEC) 


Figure C-15. Maximum contribution performance 
with only three CMGs. 





POINTING ERROR (10 4 DEG) 


Additional runs in which other CMGs were failed indicate that the 
HEAO performance specifications can be met using any three CMGs. In each 
case the CMGs must be initialized to a new null state. If. however, a full 
orbit under worst conditions must be attained before CMC desaturation or 
gyro- hang- up, the CMG momentum must be increased to 500 ft-lb-sec per 
wheel. Runs with the momentum per CMG raised to 500 ft-lb-sec illustrate 
that the gimbal angles stay relatively small over one orbital period, even with 
one CMG failed, and both pointing and jitter specifications are more than 
satisfied. If more than four CMGs are used, the induced cross coupling 
effects (see Figure C-12 at 4600 seconds) become less and the maximum 
contribution steering law performance is enhanced. The MC offers maximum 
growth potential because more CMGs can be added without altering the basic 
form ol the steering law or the mathematical manipulations. 

With the addition of magnetic torquers for CMG momentum manage- 
ment, the pointing improves by a factor of about 100 and the jitter by a factor 
ol 4. As illustrated in Figure C-17, the maximum peak values are 0.022 arc 
min pointing, 0. 16 arc min roll, and 0. 19 arc sec/sec jitter. The gimbal 
angles stay less than 4 degrees over the orbit and the magnetic system dumps 
about 400 ft-lb-sec momentum. After failing any CMG. the magnetic system 



Figuie C 1*. Maximum contribution performance with 
magnetic momentum dump. 
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will automatically force the CMGs to a new null position. No reprogramming 
or changes are required in the MC steering law. After finding a new null, the 
three-CMG system with magnetics performs the same as the four-CMG 
system shown in Figure C-17. If two CMGs are failed, direct magnetic 
torque must be used in addition to magnetic momentum dumping for control 
of a vehicle axis. Most of the time, HEAO performance requirements can be 
maintained. However, there are short time intervals of about 50 seconds 
during some orbits that 1 arc min pointing is exceeded. The MC steering law 
with magnetics is the simplest way of providing fail operational capability. 

In run MC-102, CMG number 3 and number 4 were failed by setting 
their momentum to zero. No changes were made in the steering law, nor 
were the remaining two CMGs set to a new null position. Without magnetics 
the program diverged. There was an initial momentum of 250 ft-lb-sec on 
both the Y and Z axes that produced initial vehicle rotations. The two opera- 
tional CMGs could not correct the situation without additional torque from 
either the RCS or bar torquers. In the next run (MC-103) , the magnetic 
system was used to provide both direct torque and CMG momentum manage- 
ment. The performance is shown in Figure C-18. During the first quarter 
orbit, the CMGs are automatically driven to a new null position by the magnetic 
system, during which time the initial pointing error is about 0. 16 degree and 
the roll error peaks at 15. 6 arc min at 1000 seconds. The new null position 
is with CMG angles aq - -59 and q 2 = 59 degrees. At these angles the 
CMG momentum vectors are opposite each other on the intersection of their 
momentum planes. The axis that cannot be controlled by the two CMGs is 
approximately aligned with the momentum vectors at their null position. The 
magnetic system must provide torque about the uncontrolled axes. After the 
new null is attained, the peak errors are 0.48 arc min pointing and 0.36 arc 
min roll. The jitter is well within the bounds specified for HEAO. The CMG 
gimbal angles are shown in Figure C-19. Dashed lines represent the null 
position for the two operational CMGs. The angular excursions from the 
null are less than 10 degrees once it has been established. The accumulated 
CMG momentum oscillates proportional to the gimbal angle deviations with 
peak values of about 100 ft-lb-sec (not shown). 

The maximum contribution is the only steering law that provided fail 
operational capability when any two CMGs w'ere failed. No logic is required 
to detect the failures and no modifications to the steering law are necessary 
with magnetics. The magnetic system permits any CMG-out type failure 
without any detection and logic required to define a new CMG null position. 

With the MC steering law, each gimbal is commanded independently of the 
others; therefore, no modifications are required when any CMG fails. 
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figure C lb. Maximum contribution performance with 
two CMGs tailed and magnetics. 



figure C 19 . CMC gimbal angles for the MC law with 
two CMGs failed and magnetics. 
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Pseudo Inverse Steering Law 


Typical performance with all CMGs operating using the pseudo inverse 
(PI) steering law is shown in Figure C-20. The maximum deviations are 
0 . 062 arc min on both the experiment and solar pointing axes and less than 
0.004 arc sec/sec jitter which clearly meet HEAO-C specifications. No 
momentum has been dumped. As a result of secular g. g. momentum, the 
CMG gimbal angles get rather large. As shown in Figure C-21, gimbal 
angle number 3 attains a value of 142 degrees at the end of one orbit, T^ , 

while CMG gimbals one and two reach a magnitude of 82 degrees. Although 
not shown, saturation is attained near two orbits with 0 ^ = -90, q 2 = 180, 
a 3 = 90, and 04 = 0 at which time control is lost. 

Using the pseudo inverse steering law, the vehicle angular rates stay 
very small, consequently the pointing performance is very smooth. The 
pointing error (Fig. C-20) has the same shape as the gravity gradient dis- 
turbance torque (Fig. C-3). The momentum accumulated by the CMGs is the 
same as the gravity gradient momentum shown in Figure C-4 but opposite in 
sign, and total magnitudes are identical until the time of CMG saturation after 
which the two diverge and pointing control is lost. CMG saturation also 
corresponds to gyro hang-up and to a mathematical singularity in the steering 
law algorithm. Using the pseudo inverse, those gimbal angle combinations 
which produce gyro hang-up will also cause the determinant of the ICC* I 
matrix to go to zero. Although there are an infinite number of CMG gimbal 
conditions that can produce gyro hang- up internal to the maximum momentum 
envelope, these conditions were not encountered under normal pointing condi- 
tions. But by replacing the cyclic g. g. torque by a properly directed constant 
torque, gyro hang-up could always be encountered. Much more work is 
needed to characterize the gyro hang- up problem associated with SG CMGs 
and to assess its impact on vehicle pointing performance. 

With the addition of electromagnets for continuous momentum dump, 

the performance improves by about a factor of four, as shown in Figure C-22, 

With the magnetic torquers sized at 0.4 ft— lb/ gauss and a magnetic loop gain 

K equal to 0.01 sec -1 , the maximum pointing error is 0.019 arc min. 
m 

Since momentum is being continuously dumped, the stored CMG 
momentum stays near zero. Consequently, the CMG gimbal angles, shown 
in Figure C-23, deviate less than 4 degrees from their null position. Both 
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CMG GIMBAL ANGLES (DEG) POINTING ERROR 6(10 ^DEG) 


RUN PI -74 



1 2 3 4 5 T 0 

ORBIT TIME (10 3 SEC) 


Figure C-20, Pseudo inverse pointing performance. 
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the gravity gradient and the accumulated CMG momenta are shown in Figure 

C-24. The difference between the curves represents the momentum that has 

been dumped by the magnetic system, about 400 tt-lb-see per orbit. The coil 

dipoles, shown in Figure C-2f>, have been hard limited lo 0„4 ft-lb,/ gauss and 

commanded proportionally to the vector components of the stored CMG 

momentum. The dipole commands were derived in Appendix A. As noted. 

the X -axis component saturates during the orbit. By decreasing K to 
v m 

0.001. the coils do not saturate, but not as much momentum is dumped and 
a corresponding increase occurs in CMG gimbal angles. The pointing per- 
formance is relatively unchanged. However, by increasing K to 0.1, all 

tht i colls reach saturation values. Both the stored momenta and gimbal 
angles stay near zero, but the pointing performance is degraded slightly. 

Using still higher K values causes the coil dipoles to react in a bang-bang 

manner that produces a magnetic torque which greatly degrades the pointing 

performance. Based on several runs in which K was varied, a magnetic 

m 

loop gain of 0.01 is recommended for continuous momentum dump. Figure 
C-26 illustrates the magnetic torque applied to the spacecraft as the result of 
the coil dipoles shown in Figure C--25. The components of the magnetic torque 
are about equal to that of gravity gradient in both magnitude and shape, 
indicating that the magnetic system is, indeed, counteracting the environ- 
mental torque, leaving the CMGs with relatively little to do. 

’ iir ; <-i; . n.ji m-o- ul lJu i pc ri ; .'uuu axes i long vehicle axis) is shown 
m ; m - re i. T<*. 'A iliiout magnetics, the imii error peaks at about 0.072 arc 
min. bat aitii magnetics it increases to 0.1 arc min. These data are in 
coni- a a i»; those gene rally observed when comparing performance with and 
without the magnetic system. In all cases where errors about two axes were 
root-sum-squaroii Lo get the pointing 1 error, tin* use of continuous electro- 
magnetic momentum dumping improved performance. 

CMG number 2 v\as failed by setting its momentum to zero. However, 
no change was made in computing the gimbal rate commands for the remaining 
CMGs. and the elements of the column vector in the CMG torque matrix 
corresponding to the failed CMG were 4 not set to zero. In run lT-Ga the gimbal 
angles \vr;c not initialized to a new null position 1 position at which the CMG 
momentum is zero) tor the flneo operational CMGs, The effect was to pro- 
duce a biased momentum component of 250 ft -lb-sec on the positive Y-axis. 
Without magnetics, pointing control was lost after 1000 sec orbit time and 
did not recover. The CMGs had to be initialized to a new null position using 
the RCS system alter which control was maintained until CMG saturation. 
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ROLL ERROR (10‘ 3 DEG) MAGNETIC TORQUE (FT-LB) 



Figure C-26. Magnetic torque components dumping 
CMG momentum. 



Figure C-27 0 Roll error with and without magnetics. 
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In an attempt to find a new null with CMG number 3 failed, maximum 
use was made of the magnetic system by setting = 1 for direct magnetic 

torque control and K =0.01 for CMG momentum dump. The results are 

m 

summarized by Figures C-28 through C-31. The magnetic system forced the 
CMGs to a new null position while maintaining HEAO-C pointing requirements. 
The peak pointing error (Fig. C-28) was 0. 084 arc min and the jitter was 
0. 72 arc sec/sec. The new null position was attained after about 2000 seconds 
(Fig. C-29) with <y 1 = 0, o 2 = 56.4, and -56. 4 degrees. As shown in 

Figure C-30, the initial CMG momentum was 250 ft-lb-sec but the magnetic 
system rapidly reduces it to less than 40 ft-lb-sec after 2000 seconds time. 

In so doing, the magnetic system was continuously exercised and all the dipoles 
were saturated at 1000 seconds (Fig. C-31), but toward the end of the orbit 
were operating in their linear regions. In addition to dumping the initial 
momentum of 250 ft-lb-sec, the accumulated g. g. momentum of about 400 
ft-lb-sec has also been dumped. During this orbit, the magnetic torque com- 
ponents which were generated attained magnitudes of 0. 22 ft- lb, a value 
greater than the gravity disturbance torque. With the magnetic system , the 
CMGs continued to operate about the new null with less than 4 degrees gimbal 
variation over the next several orbits (not shown) . The pointing performance 
was only slightly degraded from that shown in Figure C-22 with all four CMGs 
operating. 

In additional cases (runs PI-98 and -99) the elements of the torque 
matrix corresponding to the failed CMG were set to zero. The performance 
of the three operational CMGs equaled that obtained with all four operating 
once the new null had been established. With- one CMG out, there are three 
remaining gimbal angles and three equations that relate them to the commanded 
torque, in which case an exact inverse can be used to obtain the CMG steering 
law. Additional cases (not shown) were programmed with exact inverses for 
three CMG configurations. Data from these cases were compared with those 
of the pseudo inverse with one CMG out. The results were identical. That is, 
with one CMG out, the pseudo inverse steering law reduces to an exact inverse. 
Again, the simulation data agree with the theory of the pseudo inverse [8 and 
12 ]. 


With continuous momentum dump, the CMGs do not accumulate much 
momentum. In a subsequent run (not shown) four 25 ft-lb-sec CMGs were 
used to maintain satisfactory performance. In general, the electromagnet 
aligned with the axis of minimum inertia does more work than the dipoles along 
the other axes. In the next run (not shown) the X^-axis coil was failed, as 

well as CMG number 3, and the pointing performance was relatively unchanged. 
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Figure C-28. Performance while finding a new CMG null 
with CMG number 3 out. 



Figure C-29. Magnetics force CMGs to a new null 
with CMG number ;3 out. 
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COIL DIPOLES (FT-LB/GAUSS) CMG MOMENTUM (FT-LB-SEC) 


/ 

/ 



ORBIT TIME <10 3 SEC> 

Figure C-30. Magnetics force the CMG momentum to zero 
as a new null is found. 



Figure C-31. Coil dipoles with CMG number 3 failed 
versus orbit time (10 3 sec). 
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However, only about two-thirds of the accumulated gravity momentum was 
dumped. The system could saturate in about four orbits if a worst-case 
attitude hold were maintained with one CMG and one coil failure. 

With the magnetic system dumping momentum, two CMGs were failed 
and, without modifying the steering law, control was lost. With two CMGs 
failed, the remaining two CMGs can provide only two-axis control. The third 
axis must be controlled by magnetics or the RCS. Moreover, the pseudo 
inverse must be reprogrammed. Since there are only two unknowns (gimbal 
rates) and three known quantities (commanded torque components) , the 
pseudo inverse has the form 

c + (OCJ-'C* 
and the gimbal rates are 


a 


C~ T 

c 


In this ease, the torque matrix C is a 8 by 2 matrix and the pseudo inverse, 

C , is 2 by d matrix. To obtain maximum use of the magnetic system, both 

diiect toi quo (K 1* 0) and momentum dump (K — 0. 01) commands were 
c m 

used to drive the coils. The pointing performance is shown in Figure C-M2 
with CMC'*: number 2 and number failed. The pointing error peaks at 2350 
seconds with a value oi 0.22 arc min. In this particular ease, the roll error 
about the experiment- axis is 4.8 are min. The magnetic system is providing 
loll control and at 2750 seconds the earth’s tield is unlavorable for roll con- 
trol (Fig. C-7) . Both the Y and Z axes components (B and B ) are near 

v v y z 

zero and the axis component (13 ) cannot be used to produce roll control 

toique. Therefore, at 2750 seconds, the roll axis is not being controlled and 
the roll error is building up. However, this condition lasts only for about 
50 seconds during the orbit before the earth’s magnetic field changes and roll 
control is re-established. Although the roll error is large compared to the 
pointing error, it is still within the 5 arc min goal baselined for HEAO-C. 
Additional runs were made with only two CMGs operational and augmented with 
magnetic torques. In most cases, the HEAO-C pointing requirements were 
maintained. In those cases where the errors exceeded requirements, the 
excess errors were only for very short time intervals during the orbit, less 
than 100 seconds per orbit. It is concluded that, in all but the most unfavorable 
pointing orientations and orbit conditions, two CMGs augmented with direct 
magnetic torque can meet the HEAO-C pointing requirements. 
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Figure C-32. Performance with CMG number 2 and number 4 out 
using the pseudo inverse with magnetic control. 

With CMG number 2 and number 4 out, the momentum vectors of CMG 
number 1 and number 3 are equal but opposite in direction. The two- CMG 
system is, therefore, at a null position. The variation in gimbal angles, 
shown in Figure C-33, is less than 10 degrees over the orbital period and the 
accumulated momentum (not shown) is less than 100 ft-lb-sec. 


The Bendix Three-Gimbal Inverse 


The three-gimbal inverse steering law proposed by Bendix [6] is more 
complex than any of the other laws. Consequently, more time and effort were 
required for simulation. The four 3 by 3 matrices were inverted by several 
methods: (l) a subroutine for matrix inversion, (2) direct programming for 
each inverse, and (3) an iterative technique for which the matrix elements 
were updated at each time step. Each method gave somewhat different results. 
Ail, however, produced similar trends without magnetic dump. At about 1200 
seconds in the orbit, the torque vectors for CMG number 1 and number 3 
became colinear causing the determinants of matrices A 2 and A t to approach 
zero. Even though a three-dimensional space was still spanned by CMG torque 
vectors number 2 and number 4 and a vector aligned with number 1 and number 
3, the program diverged due to mathematical singularities. For the three- 
gimbal inverse to perform properly, it is absolutely essential that singularity 


112 



6 r 



Figure C-33. CMG gimbal angles with CMC number 2 and number 4 out 

and magnetic control. 

detection and avoidance schemes be incorporated into the steering law. How- 
ever, due to the complexity of such schemes, this was not done during the 
sim .nation. 


With the directly programmed inv erse, the singularities occurred at 
1200 seconds. In contrast, the subroutine for matrix inversion used an itera- 
tive method and, hence, was not quite as sensitive to singularities. With this 
subroutine, the program diverged at 1700 seconds. However, with the itera- 
tive technique utilizing matrix element update, the determinants were not used 
in obtaining the inverse. With this inverse, oscillations occurred at 1800 
seconds after which control was regained and pointing performance maintained 
over most of the orbit. Figure C-34 shows pointing performance using the 
iterative technique. The peak error of 5.4 arc min occurs at a time (1900 
seconds) just after the time the matrices would have singularity. As illus- 
t rat 3 d, the performance is not smooth and, many times during the three- 
fourths orbit period, the jitter exceeded HEAO requirements. 

The gimbal angles are shown in Figure C-35. Notice the sharp breaks 
where the singularities occurred (1900 seconds). The movements are irregu- 
lar especially at 1700 and 4000 seconds, although the general trends are simi- 
lar to those obtained with the other steering laws. 
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With electromagnets inserted and used for continuous CMC momentum 
management, the gimbal angles stay small, thus, avoiding the singular condi- 
tion. In these cases, the performance of the three-gimbal inverse was com- 
parable with that obtained bv the other steering laws and, therefore, is not 
shown. 


The shortcoming of the three-gimbal inverse is that most of its singu- 
larities are self- induced. That is, the basic law itself permits mathematical 
singularities that are not singularities for the other steering laws. However, 
with a suitable singularity detection and avoidance scheme, the full momentum 
envelope could probably be utilized for control purposes. Nevertheless, the 
complexity ol the scheme with the associated matrix inversion procedures 
appears to prohibit its use. 


Transpose with Torque Feedback 

The transpose type steering law is derived by taking the transpose of 
the CMG torque matrix as an approximation for its inverse. The gimbal rate 
commands are: 


r v - 


-(c T + C T + C T )/H 
11 cx 21 c.v ;;i ez 


a = 


< C 12 T cx + C 22 T c.v + VcJ /U 


rv 


-(C T : C T - C,.,T- )/H 

!•» cx 2.i c \ >iij ez 


and 


, - -(c t + c r r + c t ) /h 

4 14 cx 24 cy 114 ez' 


whcic C_ ai c the elements oi the torque matrix, T^ ^ are the components 

ol the commanded torque, and H is the momentum per CMG. Typical perl'or- 
mance is shown in Figure C-:i(i. The pointing error is u. b arc min, and the 
loll eritn exceeds that specilied lor HHAO. However, by doubling the feed- 
back gain on the roll axis (axis of minimum inertia) , the roll error (not shown) 
was less than the pointing error. The peak pointing errors are caused by two 
f.ictois: (l) Just before 4500 seconds, the CMG gimbal positions (Fig. C-37) 
were ver, near the gyro hang-up position of a { = -90, <v 2 = 0, = 90, and 

- 0 aim (2) The nonlinear terms in the Euler equations added significantly 
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Figure C-36. Performance of the transpose steering law. 



Figure C-37. CMG gimbal angles versus orbit time (10^ sec) 
for the transpose steering law. 

to the error. These terms are products of the bod} 7 angular rate multiplied by 
the accumulated CMG momentum. The body rates (jitter) peaked at about 0. 98 
arc sec/sec with over 400 ft-lb-sec accumulated in the CMGs which produced 
errors through their vector cross product terms in the Euler equations. 
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ROLL ERROR {1(T Z DEG) 


A magnetic system was added for CMG momentum dump. The perform- 
ance (Fig. C-38) was about 0.026 arc min pointing and 0.2 arc min roll. The 
peak jitter (not shown) was 0.27 arc sec/ sec. This performance is not quite 
as good as that obtained with the other steering laws using magnetics, basically 
because the torque produced per unit torque commanded is not unity. The 
gimbal angles stay less than 4 degrees when the magnetic system is used for 
CMG momentum management. During most runs simulating normal pointing- 
modes, the CMG gimbal rate limits are never attained. 



2 3 4 5 T 0 

ORBIT TIME (10 3 SEC) 


Figure Transpose* performance with magnetics. 

The transpose with torque feedback steering law is obtained by feeding 
back the CMG torque at each instant of time and subtracting it from the torque 
commanded. The commanded torque in only the gimbal rate equations (steer- 
ing law) are replaced bv 

1 ‘ i - h 

c/\ ex X 

T T = T - h 

cy cy y 

and 

T * -■ l - h 

ez cz z 
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where h is the CMG torque vector components. Figure C-39 shows the 

X y y > z 

pointing performance. The response is oscillatory during the first half-orbit 
between 0. 005 and 0. 018 degree. The peak error of 4. 7 arc min occurs at 
4700 seconds, at which time the roll error is 1. 05 degrees and jitter is 0.47 
deg/ sec. These data are the worst obtained with any steering law. On exami- 
nation of the gimbal rates they are all chattering between their hard limits of 
1 deg/ sec. Although it appears that the oscillations have a period of 100 
seconds, that could be a false conclusion since the data were hand-plotted from 
computer printout at each 50 seconds of the orbit. Therefore the gimbal angu- 
lar rates and pointing errors could be chattering at a higher frequency between 
the 1000 and 3400 second period. The CMG torque feedback has the effect of 
greatly increasing the control system gain through the CMG loop, by adding a 
lead to the system. To properly stabilize the system, tachometer dynamics 
or a lag filter needs to be inserted in the CMG torque feedback loop. However, 
this was not done during the simulation, so the results obtained are not repre- 
sentative of a properly operating system. The gimbal angles (not shown) 
approximately equaled those shown in Figure C-37, but were not as smooth. 



Figure C-39. Transpose with torque feedback performance. 
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Magnetics were added on the next run. As in the other magnetic runs, 
the gimbal angles stayed small, less than 4 degrees, over the orbital period. 
However, the gimbal rates still chattered continuously but the performance 
was greatly improved. Peak errors were 0.16 arc min in pointing, l.G arc 
iuin in roll, and 9.3 arc sec/sec in jitter. Data from this run are not shown. 


The BECO H-Distribution 

Typical performance using the BFCO steering law is shown in Figure 
C p). The peak pointing error is 0. 10 arc* min at 4800 seconds and the peak 
mil error is 0. 05 arc* min at 1000 seconds. Compared with the pseudo 
inverse* < Figs. C-'M and C-iil) , the performance is degraded by a factor of 
to. However, in additional runs (not shown) . the gain iactor on the fourth 
gimbal rate was increased from 0.0000 in run B-10 to 0.001 and the per- 
formance was approximately that obtained with the pseudo inverse. In run 
B-10. the gimbal angles t Fig. C-41) attain slightly larger values than those 
shown in Figure C-21. At the end of one orbit. O j -100. e>> - 110. a - 166, 

and a 4 18 degrees. As the run continued into the second orbit, a 

singularity was encountered at 8000 seconds, about one and one-hall orbits, 
at which time control was lost. For this run, pointing control was lost about 
one-third orbit sooner with the BFCO law than with the pseudo inverse. 

To test the ability of the BFCO steering law for distributing momentum, 
,i; . w wiviiivi-il f ro in the simulation ;iml a constant 

:; # a it-Lb was placwi on the V -axis of the vehicle. Since there is a 

! i i o < 1 1 ( * 1 1 1 ; t m i p« a cni jai oi sue i !. - 1 1 >- s > ( ‘ i»u c.u'li . ixis, (CMC saturation should 

occur at HUH.) seconds if gyro hang- -up is avoided. In this particular situation, 
tho gyro hang-up condition lu bo avoided is expected to occur at 1 oOO seconds 
with n ; on. : i a r an. and « o degrees, ihis CMC orientation 

pliKH. .‘j a. [ i the CMC.- I (;i < jiu. vec tom ou me \ / plane. Since the disturbance 

torque is pt a pi ■ 1 1 d 1 1 * *a 1 .■ • i ■ 1 1 > lias plane. the. 1 e om m a nded l ( ) r q u e will 1 ( )1 J o\v the 
disturbance so that the vector scalar product between the commanded and 
CMC torque will be zero, hence gy ro hang-np. Figure C- i 2 illustrates the 
gv'rc.) hmig-up condition ml at Ha- to Idm applied constant torque. With a constant 
torque input ( J ig. t • jol, ic ; ptuaiinu error is also constant at o.nsa are min 
until i Md saturation or gyro r m;, no. 

With the* BFCO law. saturation occurred at about 4000 seconds, as 
predicted, after which control was ioM. r l he gyro hang up condition was 
avoided by di stributing liv aluminium hmweai axes, lienee, the BFCt > 



CMG GIMBAL ANGLES (DEG) POINTING ERROR, p <10‘ J DEG) 



Figure C-40. Performance with the BECO steering law. 



Figure C-41. CMG gimbal angles with the BECO steering law’. 


120 




H-distribution law performed as advocated. The same conditions were run 
with the pseudo inverse law but, as shown in Figure C-43, the gyro hang-up 
condition was encountered at about 1500 seconds and control was lost. The 
system did not recover but stayed in the hang-up condition. The corresponding 
gimbal angles are shown in Figure C-44 for the pseudo inverse and in Figure 
C-45 for the BECO steering laws. The BECO law forced the fourth gimbal 
angle to move 180 degrees so that the full momentum envelope was utilized. 
Prior to the gyro hang-up condition, the gimbal angles were the same for 
either case. However, at 1300 seconds, the BECO law forces gimbals 
number 4 and number 2 from their null positions to avoid the hang-up condi- 
tion. In so doing, gimbal number 1 is rapidly forced past the critical 90- 
degree point, attains a peak of slightly more than 140 degrees, and then 
decreases to 90 degrees at saturation. The CMGs are saturated in the 
Y^-axis at 4000 seconds with « ] - 90, CV 2 = 0, a 3 = -90, and o 4 = 180 degrees. 

As previously defined, saturation represents the ultimate in gyro hang-up but 
cannot be avoided by any steering law, unless momentum is dumped. Control 
is always lost at saturation but some cases of gyro hang-up internal to the 
momentum envelope can be avoided by the steering law. The BECO law may 
be directional and needs further development to prove its ability in preventing 
gyro hang-up conditions. 


SUMMARY 


During the study, several CMC steering laws were evaluated. As long 
as the gimbal angles stay less than 90 degrees, almost any steering law can 
meet the HEAO requirements. With 250 ft-lb-sec CMGs, the gimbal angles 
get large within an orbit, thereby, ruling out the use of a constant gain steer- 
ing law. When one CMG has failed, the remaining CMGs must work harder. 
With several of the steering laws (the Bendix three-gimbal inverse, for 
example) the failure must be identified and corrective changes made. After 
making any required changes, the resulting steering law must be identical 
to the exact inverse of the 3 by 3 torque matrix to prevent degradation in 
pointing performance. Based on both simulation results and mathematical 
theory, the pseudo inverse steering law reduces to an exact inverse when any 
CMG is arbitrarily deactivated. With the pseudo inverse, system performance 
is not degraded by using only three CMGs for control. 
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CONTROL LOST AFTER 1500 SEC 



Figure C-44. Gyro hang-up at 1500 sec with pseudo inverse 
steering law and constant torque. 



Figure C-45. Saturation at 4000 so*:* with BKCO 
steering law and constant torque. 






To allow more than one CMG failure without affecting performance, 
consideration should be given to using more CMGs but with each sized to a 
lower momentum capacity. For the same total momentum capacity, six 
125 ft-lb-sec CMGs would permit three failures without degrading performance. 
However, with three failures, the momentum would have to be dumped about 
each half-orbit to prevent saturation under worst-case environmental 
conditions. 

Based on estimated impulse requirements lor CMG momentum manage- 
ment, the fuel weights for an all-RCS dump could become prohibitive, espe- 
cially for a growth version of HEAO, for w'hich the inertia distribution becomes 
less favorable. The alternative system recommended for HEAO utilizes 
electromagnetic torquers reacting with the earth's magnetic field to dump 
accumulated momentum. In this case, an RCS is not needed after the OAS 
burn control period. Tradeoffs show that the magnetic system is bettei than 
RCS from both a weight and reliability viewpoint. Moreover, the low torque 
levels of a magnetic system permits continuous momentum dumping without 
interfering with experiment pointing. Simulations show and analysis has 
proven that, as an added bonus, a magnetic CMG desaturation system improves 
pointing performance by providing integral control of the attitude error signal 
through the magnetic loop. Since CMG momentum is continuously dumped, 
the gimbal angles stay small (less than 4 degrees for four 250 ft-lb-sec CMGs) 
and, typically, the stored momentum is less than 20 ft-lb-sec. Hence, with 
a magnetic CMG desaturation system, the CMG momentum per wheel could 
be reduced considerably as compared with the present baseline size. Altei- 
natively, a greater depth of CMG failures could be tolerated without degrading 
performance. 


Magnetic momentum dump always keeps the CMG gimbal angles and 
momentum small, permitting linear operation of the steering law. The per- 
formance of any steering law is enhanced by the magnetic system. With 
small gimbal angles, the performance obtained by various steering laws was 
comparable. However, once the gimbal angles get large, the performance is 
usually degraded by cross coupling and nonlinear effects in the Euler equa- 
tions. Only the pseudo inverse and H-distribution laws performed without 
degradation with large gimbal angles. Moreover, maneuvers were commanded 
with the CMGs near a saturation condition to illustrate the transfer of momen- 
tum from one spacecraft axis to another — one ol the reasons lor selecting a 
near spherical CMG momentum protile tor HEAO-C. 

With continuous magnetic dumping, the constant gain steering law 
meets all HEAO pointing requirements and would be the simplest to implement. 
The maximum contribution with magnetics offers fail operational capability. 
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Even with two CMG failures, no program changes are needed. The pseudo 
inverse must be reprogrammed for two CMG failures. Neither the three- 
gimbal inverse nor the H-distribution steering laws permit two failures. 
Although the transpose with torque feedback is similar to the MC, CMG 
failures were not simulated and additional studies are needed to obtain the 
proper stabilization networks for optimum performance. For the greatest 
depth of failures without any program modifications, the maximum contribu- 
tion steering law could be used for HEAO with a magnetic system utilized for 
continuous momentum management. 


Without magnetics and requiring at least one orbit of CMG control 
prior to RCS dumping, a sufficiently large momentum envelope must be 
available for control during attitude hold modes. Over extended periods 
between dumps, the CMG gimbal angles and stored momentum become large, 
hence, cross coupling and nonlinear terms in the Euler equations can produce 
significant pointing errors. Only the pseudo inverse and H-distribution 
steering laws permitted full utilization of the CMG potential without unde- 
sirable side effects. Moreover, either law provides growth potential for 
greater pointing accuracies than is required for HEAO. But since the pseudo 
inverse also provides fail operational capability for one CMG out, without 
reprogramming, it is recommended for use on HEAO with a RCS used for 
periodic momentum management. 

In gen ci al with disturbance torques acting on all three axes, gyro 
hang-up was not encountered using pseudo inverse steering. However, an 
increase in pointing errors was observed whenever the gimbal angles were 
near a gyro hang-up or singular condition. With the pseudo inverse steering- 
law as implemented, gyro hang-up also corresponds to singularity. An 
alternate implementation ol the pseudo inverse steering law is possible that 
will completely remove singularities although internal hang-up conditions can 
still be encountered. None of the laws simulated were designed to avoid hang- 
up. However, the BECO law looks promising, although additional work is 
required to prove its worth. As used in this report, singularities are a 
mathematical occurrence which is inherent to a specific steering law formu- 
lation. Whereas, gyro hang-up is a physical orientation of the gimbal 
positions which prevents the desired torque from being produced. Currently, 
it appeals that only about 50 to GO percent of the total momentum envelope of 
single gimbal CMGs is usable before encountering a possible hang-up 
position. Much more research is needed to understand and devise ways of 
avoiding gyro hang-up conditions and is outside the scope of this report. 
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